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Abstract 

We investigate in a covariant manner the spin-induced non-geodesic motion of massive spin-| 
particles in an arbitrary gravitational field for trajectories that are initially geodesic when spin is 
ignored. Using the WKB approximation for the wave function in curved spacetime, we compute 
the 0{h) correction to the Wigner rotation of the spin-i particle, whose 0(1) contribution is zero 
on timelike geodesies. We develop conditions for the motion of observers in which the Wigner 
rotation is null. For the spherically symmetric Schwarzschild metric, we consider specific examples 
of particle motion in the equatorial plane for (i) circular orbits and (ii) radially infalling trajectories. 
For the former case we consider the entanglement for a perfectly anti-correlated EPR entangled 
pair of spins as the separate qubits traverse the circular orbit in same direction. 
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I. INTRODUCTION 



Entanglement is an important resource for many applications in quantum information 
science (QIS) including teleportation, quantum computation, and quantum communication. 
Recently, there has been a growing interest in understanding entanglement in quantum in- 
formation science beyond the confines of its non-relativistic quantum mechanical origins. A 
excellent recent review can be found in Peres and Terno [1] and references therein. Initial 
studies concerned the behavior of quantum states, both single particle and bipartite entan- 
gled states under the action of Lorentz transformations (LT) [2] , which transforms between 
different inertial (constant velocity, zero acceleration) observers. When quantum mechan- 
ics is merged with special relativity (SR) in the form of quantum field theory (QFT), the 
state of a particle is labeled not only by its spin (or helicity) (as in non-relativistic quan- 
tum mechanics (NRQM)), but also by its momentum. These two quantities represent the 
Casimir invariants of operators which commute with the Poincare transformations (LTs plus 
translations) which underly the symmetries of flat (Minkowski) spacetime, where SR applies. 

Peres et al [3] were among the first to point out that even for a single particle, this 
could lead to an observer dependent change of the (von Neumann) entropy of the reduced 
spin density matrix when the momentum is traced out after the action of a LT. Alsing and 
Milburn [4] investigated the transformation of maximally entangled bipartite Bell states 
composed of pure momentum eigenstates, and showed that while initially colinear spin and 
momentum directions are transformed under a LT to non-colinear directions, the overall 
entanglement is preserved. The amount by which the spin of a massive particle is rotated is 
given by the momentum dependent Wigner rotation angle, discussed in Alsing and Milburn 
and reviewed in the main body of the text below. Subsequent papers [5] explicitly pointed 
out the implied consequence of the previous work that there would be an apparent decrease 
in the magnitude of a Bell measurement if the measurement was made along the Lorentz 
transformed momentum direction, but the original, untransformed value would be obtained 
if the measurement was made along the transformed spin direction. Gingrich, Bergou and 
Adami [6] pointed out that since the Wigner angle is momentum dependent (since in general, 
a LT changes the magnitude of the particle's momentum), a wavepacket state composed of 
an integration of single particle states over a momentum distribution, would have the spin of 
each component state transformed differently, according to the value of its momentum. They 
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investigated a wavepacket state for a bipartite state of two spin ^ particles and showed that, 
due to spin-momentum entanglement, the reduced two particle spin density matrix (formed 
by tracing out the momentum) had a Wootters' concurrence [7] (an exact measure of two 
qubit entanglement) that depended on the inertial frame from which the state was observed 
(i.e. upon the LT considered). Similar considerations for the transformation of photon states 
under LTs were also considered by several authors [4, 8, 9]. 

A step towards more general types of motion was considered by several authors who con- 
sidered entanglement for constantly accelerated observers in flat spacetime [10] . It has been 
well known that these Rindler observers measure a thermal flux of particles (Bose-Einstein 
for bosons, Fcrmi-Dirac for fcrmions) as they move through the flat spacetime Minkowski 
vacuum, at a temperature that is proportional to the observer's acceleration (the Unruh 
effect, which is the flat spacetime analogue of the Hawking effect of black hole evaporation). 
These author's investigated the fldelity of teleportation and other entanglement measures 
for maximally entangled states of both spin zero and spin | particles. 

An important step forward in the evolution of these relativistic investigations was made 
by Terashima and Ueda [11] who investigated the transformation of single particle and 
entangled states under arbitrary states of motion (acceleration), where general relativity 
(GR) applies. Essentially, the global inertial frames of SR (zero acceleration) are now 
replaced by local reference frames, tangent to the curved spacetime (CST) at the point x, 
for an arbitrary accelerating observer at the spacetime point x. This observer is described 
in terms of a tetrad, or four 4-vector axes, three of which describe the spatial axes of the 
observer's local laboratory at x and one temporal axis which governs the local rate at which 
his clock ticks (the observer's proper time). Inside the observer's local laboratory at the 
point X SR holds, which is an embodiment of Einstien's Equivalence Prinicple (EP). An 
observer makes measurements of a particle that passes through his local laboratory at x by 
projecting the particle's momentum onto the the four axes of his tetrad. Terashima and 
Ueda showed that as the particle moves inflnitesimally from x ^ x' m. the CST, the spin of 
a particle is transformed by a local Lorentz transformation (LLT), and correspondingly by a 
local Wigncr rotation of its spin. In general, a LLT transforms between different observers 
in arbitrary states of motion, all instantaneously at the CST point x (e.g. stationary, freely 
falling, circular geodesic, or arbitrary acceleration), as will be detailed in Section II. 

In this paper we show that a consequence of the above considerations is that the Wigner 
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rotation is measured to be null in the non-rotating, instantaneous rest frame of the acceler- 
ating particle, the so called Fermi- Walker frame (FWF). In any other reference frame, the 
observer would detect a non-zero Wigner rotation angle that is dependent upon his particular 
state of motion. If the particle is undergoing force free geodesic motion (zero acceleration) 
the FWF reduces to the freely falling frame (FFF) in which all four axes of the tetrad are 
parallel transported along the particles 4-velocity (the tangent to the particle's geodesic, 
which is equal to the temporal axis of the tetrad). It is a postulate of GR that the force 
free motion of particles (massive or massless) follow geodesies. This is true if the spin of the 
particle is assumed to be zero. However, for particles with spin, even classically spinning 
particles [12], the spin of the particle couples to the curvature and leads to non-geodesic 
(accelerated) motion. In this work, we consider quantum spin half particles whose orbits 
are initially geodesic, if spin is ignored, and consider the 0{h) corrections to their motion 
when the particle's motion is defined by its quantum mechanical Dirac current [13, 14] . We 
develop the 0{h) Wigner angle for such particle motion and investigate the implications of 
entanglement of Bell states in CST. A companion article [15] explores these considerations 
for photon states and the effects of the local Wigner rotation in CST on entangled photon 
states. 

This paper is organized as follows. In Section II we review freely falling frames and 
Fermi- Walker transported frames in general relativity. In Section III we review the Wigner 
rotation and the transformation of massive positive energy, single particle states in flat 
spacetime, while in Section IV we generalize this to curved spacetime and review the work 
of Terashima and Ueda. In the Appendix A we provided a detailed derivation of the formula 
for the local Wigner rotation angle, not provided in the previous work [11]. In section V we 
develop the Dirac equation for spin ^ particles in curved spacetime, and in Section VI derive 
a WKB approximation to its solution. In Section VII we consider the 0{h) velocity and 
acceleration corrections to the initially circular geodesic motion of particles in the spherically 
symmetric Schwarzschild spacetime (derived in detail in Appendix B) , when spin is ignored, 
and discuss the consequences for entangled Bell states on two neighboring, infinitesimally 
close circular orbits. We also discuss the spin-momentum entanglement of wavepacket states, 
and show that in CST. the Wigner rotation is also dependent upon the initial orientation 
of the particle's spin in its reference frame. In Section VIII we extend the discussion to the 
case of radially infalling geodesic motion. In Section IX we present a brief summary, and 
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our conclusions. 



II. FREELY FALLING FRAMES AND FERMI-WALKER TRANSPORTED 
FRAMES IN GENERAL RELATIVITY 

In Netwonian mechanics an incrtial frame S, i.e. the laboratory from which an observer 
can make measurements, is defined as follows: (1) chose as the origin of 5* a free particle, 
for all time and (2) at one instant of time, chose three mutually orthogonal spatial axes 
defined by the orientations of three perpendicular gyroscopes. At later times, continue to 
define the spatial axes by the directions of the three orthogonal gyroscopes. Equivalently, 
parallel transport (with no rotation) the initial directions of the gyroscopes along the straight 
line trajectory of the free particle. The directions of the parallel transported spatial axes 
can be used to define Cartesian coordinates {x, y, z) with respect to the origin. All other 
reference frames (other laboratories) that move with constant velocity with respect to this 
inertial frame, are also inertial frames themselves. In special relativity (SR), an inertial 
frame is defined in exactly the same fashion, except now the universal time of Newtonian 
mechanics has to be abandoned. That is, a particular inertial frame S defines four Cartesian 
coordinates (t,x,y,z), while a different inertial frame S', travelling with constant velocity 
with respect to S, defines the coordinates {t',x',y',z'), which are related to the coordinates 
of S by the usual Lorentz transformations. 

A. Reference frames in general relativity 

In general relativity (GR) the global inertial frames of fiat (Minkowski) spacetime have 
to be abandoned for a description in terms of local inertial frames (LIF), which arc the 
SR inertial frames valid now for only a limited range of the coordinates, both spatial and 
temporal. This is just a statement of Einstein's equivalence principle, that an arbitrary 
spacetime is locally fiat. To generalize the global inertial frames of SR to LIFs in an arbitrary 
curved spacetime, one introduces four spacetime dependent, mutually orthogonal axes ea{x), 
where the hatted index a labels the four local axes a — (0, 1, 2, 3), such that ea(a;) •eg(a;) ~ r). 
Here, = diagonal(l, —1, —1, —1) is the flat spacetime metric of SR. In component form we 
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have 

9nAx)e/{x)e-^''{x) ^ri^f,, i/ = (0, 1, 2, 3). (1) 

In the above, gfj,u{x) is the metric of the curved spacetime with hne element rfs^ = 
ga/six) dx"dx^. We will use units in which the speed of light is set to unity, c — 1. With 
our chosen metric signature (+, — , — , — ) we can define the proper time r as dr^ — ds^. 
Free massive particles (m 7^ 0) follow timelike geodesies (ds^ > 0) in the curved spacetime, 
while massless particles (m = 0) follow null or lightlike geodesies {ds"^ = 0). In this work, 
we will use the term particle to mean massive objects (electrons, protons, etc. . . ) and refer 
explicitly to massless objects (e.g. photons) when needed. 

As a comment on notation, in a specific set of coordinates, e.g. x" — {t,r,9,4>), 
eo(x) = (^e^*{x),e^^{x),e^^{x),e~'^{x)j are the components of the timelike axis in the 
coordinate basis defined by x"' ( i.e. ea{x) such that ea{x) ■ eplx) — Qapix) ) and 
e-(x) = ^e. *(a;)e. ^(x), ^(x), "^(x) j are the components of the spacelike axis e-, where 
i = (1,2,3). We will denote by e{x) the collection of all four axes into a matrix for which 
the ath row is ea{x). 

The components e-'*(a;) are called a tetrad or vierbien (four- legs) [16]. Note that the 
relationship between the orthonormal basis ea{x) and the coordinate basis ea{x) (with com- 
ponents defined by (e^)^ = S^^ using the coordinates x"), is given by ea{x) — e^"(x) ea{x). 
In the following we will need the inverse matrix of tetrads e~^{x) with components denoted 

by e^"(a;). In addition, we will also need the transpose of this matrix e~^"^(,x), where T 
denotes transpose, with components given by e"' ^^{x). The inverse tetrads satisfy the dual 
to Eq.(l) i.e. 

g^^'{x)e^%x)eJix)^r]'' (2) 

where 77"^ and g'^'^{x) are the inverse fiat and curved spacetime metrics, respectively. Note 
that in matrix from we can write Eq.(l) as e ■ g ■ = rj and Eq.(2) as e~^^ • ■ = t], 
where in the last expression we have used r]~^ — rj. The matrix is the inverse of the 
tetrad e as can be seen in component form by 

eA^)eJ'{x) = S^\ e/(x)e/(x) = 5/. (3) 

We denote Greek indices {/x, ly, ■ . ■} as world indices in the arbitrary spacetime and (hatted) 
Latin indices {a, b, . . .} as local Lorentz indices in the observer's LIF. World indices are 
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raised and lowered with the curved spacetime metric gap{x) and Latin indices are raised and 
lowered with the local flat Minkowski metric Thus, the transpose inverse components 
e°^(a;) are related to the tetrad components e^^{x) by e"^(a;) = t]"''^ e^^ [x) gy^{x) . 

A vector V is a geometric object which can be decomposed in either the coordinate basis 
{ea{x)} or the local orthonormal basis {ea(x)}, i.e. V = V°'{ ). The 

utility of the tetrad and inverse tetrad is that the observer can obtain the local values V"'{x) 
of the components of a world vector V°'{x) by projecting the world vector onto the observer's 
four local axes 

V\x) = V%x) ej^{x), V%x) = V\x) e,%x), (4) 
= e\{x)V-{x). 

For a general tensorial object T°'^^{x), we flnd its LIF components by a similar 

projection onto the observer's local axes 

T%{x)^T-''^e^%x)e^\x)e,^x). (5) 

By using this set of orthonormal axes (basis vectors), the observer has made the metric 
of his laboratory locally flat, ea{x) ■ e^{x) — 77-^. The observer can subsequently construct 
coordinates in his laboratory such that the derivative of the metric gap{x) vanishes all along 
the the geodesic trajectory of the origin of his laboratory (Riemann normal coordinates). 
These coordinates are valid only if the observer's laboratory is sufficiently "small" spatially 
and measurements arc made over "short enough" times, otherwise, these coordinates lines 
can cross each other and thus become invalid for making observations. Such a reference 
frame, in which Qapixo) ~^ Vab a-nd d^Qap^xo) — where Xq are the coordinates of the 
origin of the observer's laboratory along its geodesic trajectory in the curved spacetime, is 
called Si, freely falling frame (FFF). 

The above FFF is the local analogue of the inertial frame of SR. For the motion of free 
particles, i.e. geodesies, it is the most "natural" frame from which the observer can make 
measurements. However, in general this is not the only way to define the observer's local 
laboratory. For example, the origin of a stationary observer's laboratory which sits at fixed 
spatial coordinates would in general experience an acceleration and possible spatial rotations. 
All local laboratories can be related to each other by spacetime dependent local Lorentz 
transformations (LLT) relating their choice of orthonormal bases i.e. e^(x) = Ag^(x) %{x). 
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The LLTs are independent of the general coordinate transformations (GCT) dx'°'{x)/dx^ 
that can be made in the curved spacetime, that relate quantities (vectors, tensors, etc. . . ) 
in the same spacetime described in the new coordinates x'°'{x). Thus, the tetrad compo- 
nents e^"(a:) transforms as a contravariant world vector (index a) under general coordinate 
transformations, and as a covariant local Lorentz vector (index a) under LLTs. The lo- 
cal metric r]^ transforms as a symmetric covariant local Lorentz tensor of rank two under 
LLTs, but as a scalar under GCTs. The reverse is true for the world metric Qapix) which 
transforms as a scalar under LLTs and a covariant tensor of rank two under GCTs. As a 
computational consequence of the freedom to make arbitrary LLTs, one can always perform 
computations in GR with the components of vectors (tensors) referred to a coordinate ba- 
sis V{x) = V^lx) ea{x), and then transform to the components V^lx) with respect to an 
orthonormal basis ea{x) by a LLT as in Eq.(4) and Eq.(5). 

If u{x) — u°'{x) ea{x) is the 4- velocity of a free particle with components u°'{x) — dx'^/dT 
in a coordinate basis, the condition for force-free or geodesic motion is 

VuM"(a;) = u^{x)Vpu''{x) = u^{x) (df^u'^ix) + T\^{x)u^{x)) = 0. (6) 

In general, Vu V'(x) = {^uV^i^)) ^a{x) is the total derivative (or directional derivative) of 
the the vector V{x) along the geodesic with tangent u, which generalizes to curved spacetime 
the concept of the derivative d/dr along a curve. V;3 is called the covariant derivative and, 
when one is concerned solely with the components of a vector, Vi3V"{x) is referred to as 
the covariant derivative of the vector K"(x), which is commonly denoted as V".^(a;). In the 
above, r\^{x) is the usual metric (Christoffel) afiine connection defined with respect to a 
coordinate basis ea{x) from V/3eQ,(x) = r\p{x) ex{x), and computable from the metric via 
the metric compatibility condition ^ugapix) — 0, 

rV(^) = I g''^ {d, gp.{x) + g^,{x) - g,.{x)) . (7) 

The Christoffel connection indicates how the tetrad twists and turns as it moves from x — > x', 
by specifying the rule for parallel transport of tensorial objects in the CST. Equation (6) 
states the directional derivative of the 4- velocity /x(x) along itself is zero, which is the 
geometrical statement that the 4- velocity is parallel transported along itself. More physically, 
the components of the 4-acceleration experienced by a particle along an arbitrary trajectory 
with a general 4- velocity w in a coordinate/orthonormal basis is given by 

a"(x) = u^{x) V^'u"(x), a"(x) = a"(x)e„'^(x). (8) 
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This is the external, non-gravitational acceleration that would have to be applied to keep 
the particle on its trajectory. Thus, Eq.(6) is a statement that the origin of the observer's 
laboratory experiences no acceleration along its geodesic trajectory To construct the FFF 

of the observer's laboratory, we define the orthonormal basis as follows: (1) we choose the 
temporal basis vector e^i^x) to be the observer's 4-velocity u{x) and (2) require that the 
remaining orthogonal spatial basis vectors ei{x) be parallel transported along the geodesic, 
i.e. 

FFF condition: e^ix) = u{x), Vu ea(x) = u'^{x)Vp ea(x) = a = (6, i, 2, 3). (9) 
B. Fermi- Walker transported frames in general relativity 

So far the issue of the non-rotation of the laboratory axes along its trajectory in a curved 
spacetime has not been addressed. For general non-geodesic motion with 4-velocity u (from 
now on dropping the indication of the spacetime dependence (x) unless explicitly needed) 
a particle will experience a 4-acceleration a. Note that the 4-velocity of a massive particle 
is normalized to unity, u ■ u = 1, which follows directly from the expression for the metric 
ds^ — dr^ — Qocpix) dx'^ dx^ and the definition of the 4-velocity asu — dx/dr. Consequently, 
by taking the covariant derivative of both sides of this normalization equation one finds that 
the 4-acceleration a — Vu u is orthogonal to the 4-velocity, a-u — 0. If we use this particle 
as the origin of an observer's laboratory frame, we require that the tetrads carried by the 
observer "not rotate." Some care has to be taken in curved spacetime when clarifying the 
precise definition of the concept of the non-rotation of the tetrad, since, as the particle 
progresses along its trajectory, its 4-velocity will change, which can be considered as a 
rotation (LLT) in the instantaneous frame defined by u and a. Non-rotation then means we 
accept this inevitable rotation, but require that any spacelike vector w orthogonal to both u 
and a undergo no additional spatial rotation. For a general vector V this last requirement 
is ensured by the transport law 

Wy,V ^{V ■u)a-{V ■a)u = n{V). (10) 

The rotation tensor ft accomplishes the required task as evidenced from the following special 
case: (i) fl{u) — a, which returns the definition of the 4-acceleration, and (ii) for a spatial 
vector 11? orthogonal to both u and a, fl{w) — 0, which reduces to the statement that w is 
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parallel propagated along the trajectory. A vector evolving according to the rule in Eq.(lO) 
is said to be Fermi- Walker (FW) transported along the the particle's trajectory with tangent 
u. A tetrad e(a;) with each axis satisfying Eq.(lO) is called a Fermi- Walker frame (FWF), 
and represents the instantaneous, non-rotating rest frame of an accelerating particle. Note 
that for geodesic motion, in which a = 0, FW transport reduces to parallel transport. In 
later sections, we will see explicit examples of parallel and FW transported motion of the 
tetrad defining the observer's local coordinate frame. 

III. WIGNER ROTATION 
A. Flat Spacetime 

In flat (Minkowski) spacetime, the positive energy, single particle state of a massive 
particle forms a spinor representation of the inhomogcncous Lorentz (Poincare) group [17]. 
These states denoted by |p, cr), are labelled by their spatial momentum p (where p^ = 
E — yjp"^ + m^)) and the component a of spin along a quantization axis in its rest frame 
(typically taken to be along the third spatial direction) . Under a Lorentz transformation A 
the one-particle state transforms under the unitrary transformation U (A) as 

C/(A)|p» = D'^l^^iKp)) |A?,0, (11) 

a' 

where j is the spin of the particle, the summation is over a' — {—j, — j + 1, . . . , j) and Ap are 
the spatial components of the Lorentz transformed 4-momentum, i.e. p' where p'^ = A^^p" . 
In this work wc will be primarily concerned with spin-^ Dirac particles (j = |). In Eq.(ll), 
D-'^i^{W{A,p)) is a (2j + 1) x (2j + 1) matrix spinor representation of the rotation group 
SU{2), and W{A,p) is called the Wigner rotation angle. The exphcit form of the Wigner 
rotation in matrix form is given by 

W{A,p)^L-\Ap)-A-L(p), (12) 

where L{p) is a standard boost taking the standard rest frame 4-momentum k = (m, 0, 0, 0) 
to an arbitrary 4-momentum p, A is an arbitrary LT taking p — > A p = Ap, and L~^{Ap) is 
an inverse standard boost taking the final 4-momentum Ap back to the particle's rest frame. 
Because of the form of the standard rest 4-momentum k, this final rest momentum k can 
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at most be an spatial rotation of the initial standard 4-momentum k, i.e. k — W(A,p) ■ k. 
The rotation group 0(3) is then said to form (Wigner's) little group for massive particles, i.e. 
the invariance group of the particle's rest 4-momentum. The explicit form of the standard 

boost is given by [17] 

^0 



m 

Ti _ P' rO _ 
^ — — ) ^ i — > 
m m 

= ^'.-(7-1)0, z,j = (1,2,3), (13) 

where 7 = p^/m = E/m = e is the particles energy per unit rest mass. Note that for the 
flat spacetime metric rjap =diag(l, — 1,— 1,— l),po =P° and pi = —p\ 



IV. CURVED SPACETIME 



In curved spacetime, essentially everything above goes through unchanged except for 
the important fact that single particle states now form a local representation of the inho- 
mogeneous Lorentz group at each spacetime point x. Thus, a single particle state is now 
represented as \p'^{x), a) i = (1, 2, 3) where p'^{x) are the spatial components of the particle's 
4-momentum p = p"'{x) ea{x) in the local orthonormal basis {ea(a;)}. Since the particle's lo- 
cal 4-momentum transforms under LLTs A{x), the single-particle state transforms unitarily 
via the local version of Eq.(ll) 

U{A{x))\p\x),a) = J] D'^liWix)) \{Ap) \x),a'), (14) 

where W{x) is the local Wigner rotation given by the local version of Eq.(12), 

W{x) = L-\Ap{x)) ■ A{x) ■ L{p{x)). (15) 

Following Tcrashima and Ueda [11], let us consider how the spin changes as we move from 
one spacetime point in curved spacetime to another along an arbitrary timelike trajectory. 
Let our particle initially be at a spacetime point with coordinates x and 4-momentum 
p{x) — mu{x). At a small proper time dr later the particle has moved along its trajectory 
with tangent u to the point with coordinates x' — x -\- u{x) dr and new 4-momentum 
p{x) -\-5p{x), illustrated in Fig.(l) [18]. Since the spin of the particle is defined locally with 
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p(x) = mu(x) 



FIG. 1: The observer's local laboratory (small box with man) at the curved spacetime point x, 
defined by the orthonormal tetrad ea(x). The three spatial axes e^(a:), i = (1,2,3) are located 
at the origin of the observer's laboratory, while eQ(a;) = Uo{,s(a;) is the temporal axis, defined 
as his 4-velocity, or the tangent to his geodesic trajectory. A particle of 4-momentum p{x) = 
mu{x) and world components p°(x) passes through the observer's local laboratory. The observer 
measures the local components p"(x) of the particle by projecting p(x) onto the four tetrad axes, 
p"'{x) = e"" ^{x) p" (x) . At a small proper time later dr, the particle has moved from x"' — > x'" = 
+ u°'{x)dT, which is measured by the observer in his local laboratory at the spacetime point 

respect to the observer's reference frame, defined by the tetrad that is carried along with 
him at the laboratory's origin, we are interested in the momentum change 5p{x) relative 
the orthogonal basis vectors, i.e. 5p{x) — Sp^{x) ea{x). For small dr we expect Sp°'{x) to be 
proportional p^{x) and to dr so we seek an expression of the local 4-momentum change in 
the form of 



From the definition of the local 4-momentum in the observer's reference frame as a pro- 




(16) 



where A"-(x) is an infinitesimal LLT, 




(17) 
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jection of the world 4-momentum onto the local axes p°'{x) — ^{x)jf{x) we have 

8p^{x) = e%(x) 5jf{x) + (5e%(x) p^{x). (18) 

Equation (18) contains two terms; the first 5p^{x) is the change of the world 4-momentum 
components p^{x) as the particles moves from x'^ x'^ in the underlying curved spacetime, 
and the second 5e"^(x) is is the change in the tetrads components (here, the inverse transpose 
components) e°- ^{x) which are used to project the world 4-momentum components onto the 
observer's local laboratory axes, p'^{x) — > p"'{x). The first term is simply given by 

5p^(x) = VuP^(x) dr = u''(x) V^jfix) dr ^ma^(x)dT, (19) 

where we have used the definition of the 4-momentum in terms of the 4-velocity p{x) — 
mu{x) and the definition of the acceleration in Eq.(8). Since p • p — {mcf — rri^ we can 
write Eq.(19) as 

5^{x) = - K(x)p,(x) -p'^(x) a,{x)] p^'ix) dr, (20) 

= —[n(p)r(x) dT=— Q^'Jx) p^ix) dr, 
m m 

where Q, is the "non-rotation" matrix (i.e ensures rotation only in the u-a plane) on the 
right hand side of the Fermi- Walker transport equation in Eq.(lO). 

For the second term in Eq.(18), we note that inverse transpose of the tetrads define a 
set of 1-forms defined by e"(x) = e"^(x) dx^ dual to the orthonormal basis vectors ea{x). 
The covariant derivative of an arbitrary set of basis vectors (orthonormal or not) is given 
by Vi,ea(x) = loJ' ^{x) ei^{x), which generahzes the Christoffel connection discussed just 
before Eq.(7), defined in terms of a coordinate ea{x) basis by V pea{x) = r\i^{x) ex{x). 
Since the 1-forms e"'{x) are dual to the orthonormal basis vectors ea{x) in the sense that 
e"(x) • e^(a;) = (5"^, the covariant derivative of the 1-forms (using V,/(5"- = 0) is given by 
'V„e"'{x) — — a;^^^(x) e"(x). Expanding this out in terms of coordinate 1-forms dx^ yields 
V^e"' n{x) — — a;^"g(x) e"^(x) from which we obtain the connection coefficients as 

^/s(^) = -ei'{x)V.e\{x) = e\{x)V.e-/{x), (21) 
where the last expression follow from utilizing f e- ^(x) e"^(a;) j = V,/5"- = 0. Therefore, 
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the change in the tetrad components 5e"^(x) is computed as follows 

5e\{x) = V^e\{x)dT, 

= u''{x)Vue^^{x)dT, 

= -u^{x)uj/-^ix)e\{x)dT, 

^ x'y{x)e\{x)dT. (22) 

where we have defined local rotation matrix 

X'y{x) = -u\x)u:^\{xl 

= e.^(a;) Vue%(a;) = -e%(a;) Vu e^^(x). (23) 

Substituting Eq.(20), Eq.(22) and Eq.(23) into Eq.(18) yields 5p^{x) in the the desired 
form of Eq.(16) with the infinitesimal LLT A"-(x) given by [11] 

A«g(x) = 1 [a\x)p^,{x) -p«(x)ag(x)] +x"s(x), (24) 

fit 

where we have used V°'{x) = e"' ^{x)V'^{x) for an arbitrary vector V{x), and V{x) ■ V{x) = 
V^{x)Vn{x) — V^{x)Va{x). As discussed above, the first term in Eq.(24) arises from the 
right hand side of the FW transport law and involves the local rotation of the observer's 
4-velocity in the u{x)-a{x) plane only. The second term x"^(x) arises when the observer 
chooses not to FW transport the tetrad that defines his local reference frame. For a vector 
w{x) perpendicular to the u{x)-a{x) plane, x"g(a;) produces a rotation in the plane per- 
pendicular to u{x) and w{x), i.c x(.t) ■ u{x) = and x(^) ' it'(a^) = 0. If g'-^{x) are the 
spatial axes of a non-FW transported tetrad, and e^(x) are the spatial axes of a second FW 
transported tetrad, then the former will precess locally relative to the latter with an angular 
velocity vector equal to w{x) (see MTW, pl74 in [16]). 

We can now calculate the components of the local Wigner transformation W°'^{x) that 
appear in Eq.(14) which determines how the spin of the particle rotates locally (|p*(a;), a) — > 
|p'*(a;), cr') = U{A{x))\p^{x), a)) as the particle traverses from x^ x'^ along its trajectory 
in curved spacetime. Performing this calculation to first order in dr using 

W\{x) = 5"^ + r ^(x) dr, (25) 

in the definition of the Wigner rotation Eq.(15), and additionally Eq.(17) for the form of an 
arbitrary LT to 0{dT), one can derive (after a lengthy calculation) the following expression 
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for the infinitesimal local Wigner rotation '&°'y{x) [11] 

d%{x) = d\{x) = = 0, 

d^.{x) = X^x) + ■ (26) 

■' ■' p'^[x) + m 

A derivation of Eq.(26) is given in Appendix A. For a particle of spin-j the rotation matrix 

d'^^^{W{x)) that appears in Eq.(14) is given to 0{dr) by 

D^J,l{W{x)) ^I + i [^23(^) Ji + ^3i(^) Ji + ^iiix) Jsl dr (27) 

where [^j, ^j] = ^^ijkJk commutation relations for SU{2) with the constant fiat 

spacetime spin-j matrices {Ji\. 

For the case of spin-|, to which we now specialize, we have {Jj = |(7j} where {uj} are 
the usual flat spacetime constant Pauli matrices. The infinitesimal unitary transformation 
of the state \p^{x), a) as the particle moves from — > x'^ is given by [11] 

U{K{x))\p\x)A) = (^/+l%(x)ciT) \p\x')A) 

- \ (^3i(^) - ^^23(^)) ^tIpX^O, i), (28a) 
U{K{x)) \p\x), i) = ^ (^3i(^) + ^%(^)) dT\p\x'), T) 

+ {l-\^2^{x)dT^ \p\x')d). (28b) 

where we have used the notation o"={|,— |} = {t;i}- 

Wc can iterate the formula for the infinitesimal local Wigner rotation to obtain the 
finite rotation between an initial and final point in spacetime {x{Ti),x{Tf)}. Breaking 
up the trajectory into N infinitesimal time steps of length tij/N where r/^j = dr — 
jl' {9ixv{x) dx^ dx'^y^'^ is the total proper time between the two events, and x'^ — x^{Ti + 

N 



k=0 



Til 

N 



= Texp 



drrMr)) 



n 



(29) 



In the last expression the time ordering operator T is required since, in general, the infinites- 
imal local Wigner rotations '&"'^{x{t)) do not commute at different locations x^{r) along the 
trajectory. 
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Some immediate observations can be made from the above formulas for the Wigner rota- 
tion matrix Eq.(26), which depends on the infinitesimal LT matrix A"g(x) Eq.(24) and the 
particle's 4-momentum. If we choose the observer's local laboratory to ride along with the 

particle, i.e. by selecting the observer's temporal tetrad vector to be equal to the particle's 
4- velocity 65(2;) = u{x), the particle is observed to be instantaneously at rest. The local 
spatial components of the particle's 4-momentum are then zero since 

p^{x) =p^{x)ei^^{x) =me^^{x)ei^^{x) =m6^^, ^ p\x)=0. 

Since the non-trivial infinitesimal portion of the Wigner rotation matrix ^\{x) Eq.(26), 
depends only upon the spatial components of the local 4-momentum p'^{x), the former reduces 
to '&\{x) = (since, for this case, X\{x) = x^ji^) from Eq.(24), for the same reason). 

From Eq.(23) we have x"g(^) — ~^"^(^) e^'^(x), which vanishes, in particular, for the 
case of geodesic motion in which the 4-acceleration a{x) — 0, and the choice of a FFF 
tetrad Eq.(9). For this tetrad choice the Wigner rotation matrix Eq.(25) reduces to the 
identity matrix W'^'^{x) = i.e. from this frame the observer detects no Wigner rotation 
of the particle's spin. Note that the non-rotation of the particle's spin depended upon the 
particular choice of the observer's reference frame, i.e. the FFF tetrad for geodesic motion 
discussed above. Any other choice of the spatial tetrad vectors would yield a non-trivial 
Wigner rotation '&\{x) — x^ji^) fo^" the geodesic motion, in which x*-.(a:) describe the spatial 
rotations of the observer's reference frame relative to the FFF. 

For the case of arbitrary motion the particle, with tangent u and a = Vu u ^ 0, does 
not move on a geodesic. However, in analogy to the above discussion, a similar reference 
frame can be found for which the observer detects no Wigner rotation. The observer again 
rides along in the instantaneous rest frame of the particle, with the selection of 65(0;) = u{x) 
implying that = and '&\{x) — X^ji^)- The specific choice of tetrad is governed by 

the requirement that it is FW transported along the particle's trajectory. This leads to the 
vanishing x\{x), the space-space portion of x"g(a^)- This can be seen a follows: let V in the 
FW transport equation Eq.(lO) be any of the four tetrad vectors e^(x). In component form, 
Eq.(lO) becomes 

Vue- ^(x) = [ef (x) u,{x)] a^{x) - [e^'^(x) a,{x)] M^(a;). 
Multiplying this expression by — e^"(x) and using the definition x"g(^) — 
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in Eq.(23) yields 



X%(a;) = -[ef (x) u.{x)] [a>^{x) e/(x)] + [ef (x) a.{x)] e/(a;)]. 

Finally, substituting the expression for the choice of the temporal tetrad as the 4- velocity 
of the particle for the instantaneous rest frame of the particle, u^{x) — e^^{x) and Uu{x) — 
e^°(a:), and using the orthonormality of the tetrad vectors Eq.(3) produces the expression 

X'iix) = -5g°K(x)e/(x)] + [e/(x)a^(x)](5/ 

= -a"(x)5^° + 5o"ag(x), intheFWF (30) 

which implies X^ji^) — (^^e non-zero time-space components x'^ii^) — ~X\{^) — '^li^) 
describe local spatial acceleration (boost) in the u-a plane), and hence a null Wigner rotation 
in this FW transported reference frame (FWF for short). Any other choice of tetrad would 
lead to the observer detecting a non-trivial Wigner rotation, relative to the FWF. 

In the following sections we will investigate the 0{h) correction to the particle's 4- velocity 
and 4- acceleration when the classical geodesic motion of a spin-i particle (defined by pre- 
scribing u independent of the particle's spin such that Vu u = 0), is replaced by the quantum 
mechanical Dirac current of the particle. To accomplish this task, we must first discuss the 
Dirac equation (DE) in curved spacetime [19], which we turn to in the next section. 

V. THE DIRAC EQUATION IN CURVED SPACETIME 

The covariant derivative Vq discussed in the previous section transforms world tensors 
into world tensors under GCTs. The tetrad and inverse tetrad components allow us to 
relate the components of a world tensor T°'p^{x) to the observer's LIF components T°'-^{x) 
by Eq.(21). Therefore, it would be desirable to introduce the concept of a derivative that 
is covariant (i.e. transforms well) under LLTs. For world tensors, the affine connection 
Eq.(7) was introduced on the spacetime to define the notion of parallel transport. For a 
world vector V°'{x) at a spacetime point P with coordinates x, one defines the parallel 
translate V^"{x x + dx) at point Q with coordinates x + dx the vector V^i^x) at P 
by V^^{x ^ x + dx) = V'^{x) - V'' ^p{x)V^{x) dx^ . The covariant derivative V pV^ix) is 
then defined as the subtraction of the vector "already at the point Q", V°'{x -\- dx), and the 
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parallel translate at Q, — > x + dx), in the limit that Q approaches P 

V^ix + dx) - V^^{x ^x + dx) 

We now wish extend above concept of a covariant derivative for world vectors to local 



VpV^ix) = jini^ ^ — -f- = dpVix) + V^f^ix) V^{x). (31) 



Lorentz vectors V^x). We do this by introducing a local spin connection uj^"'-^{x) used to 

define the parallel translate V^{x — > a; + dx) at the point Q, of a local Lorentz vector V°'{x) 

at point P, by V^^{x x + dx) = V^i^x) — uj^°'-^{x) V^{x) dx^^. In the analogy with world 

vectors, the local covariant derivative X'^\^"(a;)is defined by 

V^(x + dx) - VHx -^x + dx) . 
V.V^ix) = hm ^ ^ = d,V-{x) +uj^\{x) V\x). (32) 

One could further introduce the notation T>a = e^^{x)V^, da = e-^(a;) 9^ and a;, "^(x) = 
eg^(x) a;^"-(x) so that Eq.(32) only contains LIF (Latin) indices. However, such notation 
connotes the existence of a set of local coordinates y"'{x^). Such coordinates do exist, 
centered on the origin of the observer's laboratory (e.g. Riemann normal coordinates), but 
as discussed in the previous section, have only a limit range of applicability. We will retain 
the notation to emphasize our interest in the change of local Lorentz quantities as we 
move from point to point — > + dx"' in the curved spacetime. 

By requiring that Eq.(32) is "compatible" with Eq.(31) in the sense that we can transform 
between equations using the tetrad/inverse tetrad Eq.(4), we obtain the equation for the 
spin connection given by Eq.(21) (see Lawrie [20]). Equation (21) can be rearranged into 
the following form 

D,e,^{x) ^ 9^e/(x) + r\^{x) e,\x) - u;^\{x) (x) = 0, (33) 

which defines the total covariant derivative which transforms properly under both GCTs 
and LLTs, with a Christoffel connection for every Greek index and a spin connection for 
every Latin index, respectively (note: the covariant derivative V^j in Eq.(31) acts only on 
world (Greek) indices, while the local covariant derivative T>^ in Eq.(32) acts only on local 
(Latin) indices). Equation (33) can be considered a compatibility requirement of the tetrad 
(the first veirhein postulate, see Ortin, [21]) analogous to the metric compatibility equation 
V^5'Q,/3(a;) = which defined the world Christoffel connection. 

In general, the spin connection can be defined completely in terms of the (orthonormal) 
tetrad defining the observer's LIF, via 
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where 

[ea(x), e^ix)] = -2n.^%x) e,(x), ^^^(x) = e/(x) e^^^ a[^e%](x) (35) 

define the Ricci rotation coefficients which are the commutators of the basis vectors e^{x) = 
e^'^{x) d^. A non-holonomic frame is one in which the fls do not vanish, while a coordinate 
basis is one in which they do (i.e. e-'^(x) — S^'^). 

So far the above discussion has been in terms of vectors and tensors. The utihty of 
the tetrad formahsm is that it allows one to introduce spinor (integer and half-integer) 
representations of the Lorentz group through the spin connection (so named), which is 
necessary in order to describe fermions. This is the only known method by which to describe 
spinors in curved spacetime in arbitrary coordinates, and thus the only known method to 
couple fermions to gravity (see Ortin, [21]). If we denote a general spinorial quantity in the 
LIF of spin-j by '^^{x) with A — (— j, — j + 1, . . . ,_;) taking on 2j + 1 values, then total 
covariant derivative of 'if^{x) is given by 

D,^^{x)^d,^^{x) - ^a;/^(x)r(,)(Eg,)^^*^(x). (36) 

In Eq.(36) r(j)(S^g) is the matrix representation for spin-j of the flat spacetime generators 
Eg - of the Lorentz group. These are constant matrices that satisfy the commutation relations 

i^a h d] = -Vac ^hd - Vbd ^ac + Vad ^bc + %c ^ad' (37) 

For vectors we have the representation r(i)(E.^-)"- = ^rj^f- rj^-^f^, while for spinors (j = 1/2) 
we have 1(1/2) (E-^) = ^l[a,Tb] = \ [7a,7s]> where 7a = r/a^T^ are the usual flat spacetime 
Dirac gamma matrices satisfying the anti-commutation relations {7a, 7;,} = 2r7.g. These 
representations of the Lorentz group lead to the following formulas for the total covariant 
derivative of LIF vectors and spinors, respectively [21] 

D,V\x) = d,V\x)-\ujJ'\x)T^,^(Y.if^v\xl 

= d,V\x)+u^^\{x)v\x), (38) 

and 

D,^l^\x) = d,ij\x) - \ujJ'\x)V^,,,){^^,,)\r{x), 

= d,i^\x) - \u^,',{x) {b\i'])\r{xi 

= d,ij^{x) + {r,)\{x)nx), (39) 
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where we have defined the spinor connection r^(a;) for j — 1/2 

using the conventional definition of cr"^ = | [t",?''], and transforming the first index of 
(^cahi^) from Eq.(34) using the inverse tetrad components. Henceforth, we shall suppress 
the explicit spinor indices (unless needed) and write Eq.(39) as Dfj^ip{x) = {d^ + T^)ip[x). 
The promotion of the Dirac equation in flat spacetime (with global coordinates x°') 

(i'f"' da — m/fi) il){x) =0, flat spacetime 

to curved spacetime, using the minimal coupling prescription da ^ Da, becomes 

(i 7" Da - m/h) ii{x) = {i 7°(a;) - m/h) ip{x) 

— {i^"{x) [da + rQ,(a;)] — m/K) tp{x) — 0, curved spacetime (41) 
in world coordinates a;" , where we have defined the curved spacetime Dirac gamma matrices 

by 

r{x)^Yea"{x), (42) 

and we have used 7" Dg^ — 7" e^"(a;) — 7" (a;) D^. 

The FFF observer in curved spacetime can always construct coordinates such that the 
metric Qajiix) and the Christoffel symbols r^^^(a;) vanish along the geodesic trajectory (e.g. 
Riemann normal coordinates, FW normal coordinates, . . .). However, the observer cannot 
in general choose coordinates so that all the second derivatives of the metric vanish along the 
geodesic, unless the spacetime is fiat. This is described by the Riemann curvature tensor, 
which has its most direct definition in terms of the commutator of the covariant derivative. 
For vectors and spinors in the LIF this can be defined as 

[D,,D^]V\x)^R^^~^\u^{x))v\x), (43) 

and 

[D,,D^]^l;{x) ^'-R^,,i{u{x))a''^l;{x) = ^ i?^.,,(r(x)) <7^^(x) V'(x). (44) 

In the above, R^^^^{uj) is the Riemann curvature tensor formed from the spin connection, 
and is related to the usual Riemann curvature world tensor Rj^^cf^T) constructed from the 
Christoffel symbols by 

R,J(n^)) = e\{x) e/(x) R,J(u;(x)), (45) 
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where 

R,J{m) = d, r^^jx) - r^^jx) + t\^{x) r^jx) - t^,^{x) riJx), (46) 

and 

R,J{^{^)) = d, - UJ,a - ^.c + ^ua ^ ,c ■ (47) 

Equation (44) and Eq.(47) will prove useful in the next section. 

VI. WKB SOLUTION TO THE DIRAC EQUATION IN CURVED SPACETIME 
A. Wavefunction 

In this section we develop a WKB solution to the Dirac equation in curved spacctime 
Eq.(41) [13, 22]. We make no approximation to the strength or form of the gravitational 
field (metric), but instead only keep the lowest 0{K) quantum correction to quantities of 
physical interest. We develop the WKB expansion of ■0(x) as 

oo 

(48) 

n=0 

where the action S{x) is real. Inserting this expansion into Eq.(41), and using result that 
for functions DaS{x) = da S{x), we have, upon equating like powers of h 

[7"(x) daS{x) + m] M^) = 0, (49) 

[r{x) daS{x) + m] i/j.ix) = -^^{x) Da M^)- (50) 

Since the term in the square brackets in Eq.(49) is a matrix, the condition that this equation 
has non-trivial solutions requires its determinant to be identically zero 

det[7"(a;) daS{x) + m] = 0. (51) 

Equation (51) reduces to the classical Hamilton- Jacobi equation upon using Eq.(49) twice 

= r(x)daS(x)^f'(x)dpS(x) ^^{r(x),^^(x)}daS(x)dfsS(x)^g^^(x)daS{x)dpS(x) 
= d''S{x)daS{x). (52) 
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Defining the particle's world 4- momentum Pa{x) — mua{x) as the normal to the surface of 
constant action S{x) we have 



Pcc{x) = -daS{x), p°'{x)pa{x) = m^, 



(53) 



the later equation of which simply states the normalization of the particle's world 4- velocity 



u"'{x) Ua{x) = 1. 



(54) 



Thus, to 0{1) in h, the phase of any quantum mechanical particle in curved spacetime, 
regardless of its spin, is given by the classical result 



S{x) 



Pa{x) dx"". 



(55) 



The above form of the action was suggested by Stodolsky [22, 23] who pointed out that 
the action of a free particle is given by S{x) — m J ds. Writing the line element as ds — 
gap{x)dx"dx'^ /ds and defining Pa{x) = mga/3{x)dx^ / ds = mUa{x) reproduces Eq.(55). 

The determinant condition Eq.(51) and the resulting Hamilton- Jacobi equation Eq.(52) 
arise in the solution of the DE in flat spacetime. As such, in the observer's local FFF the 
general solution to Eq.(49) takes the flat spacetime form 

where P-\{x) and Pi^x) are scalar functions and the positive energy spin up ipQ (x) and spin 
down ipQ^^ (x) Dirac spinors [24] are given by 

/ ^ \ / ^ 



4'\x) 



E + m 
2m 



1/2 





p3 



E + 
\ E + m ) 



4'\x) 



E + m 
2m 



1/2 



1 • 2 
P — tp 

E + m 

3 

p 

\~ E + m J 



(57) 



with 



(58) 



{E{x), p\x)) =p\x) =p''{x)e^\x). 

Equation (58) states that, in general (i.e. not necessarily a FFF), an observer carrying local 
tetrad axes ea measures the world 4-momentum p°'{x) of a particle (massive or massless) 
crossing his laboratory by projecting p"(x) onto his local axes (in the above using components 
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of the inverse tetrad). By Eq.(l) and Eq.(2), the observer's axes form an orthonormal basis, 
so that the metric is locally flat 77-^ for the observer. Thus, the first equality in Eq.(58) states 
that the components of p"(a;) measured by the observer take the usual special relativistic 
fiat spacetime form p'^ = E = m/{l — v'^{x)) and = mv\x)/{l — v'^{x)) were v'^{x) are the 
spatial components of the locally measured velocity. 

For the particular case of the freely falling frame that we are considering, p°'{x) is the 
world 4-momentum of the observer's origin, where the spinor resides. This is the curved 
spacetime generalization of the particle's rest frame. As such, by the orthonormality of 
tetrad Eq.(3) we have 



p^{x) = {E{x), p\x)) = (m, 0, 0, 0), in a FFF. 



(59) 



As the observer traverses his geodesic trajectory in the curved spacetime, the components of 
his 4-momentum at each spacetime point will have the form of Eq.(59) in the FFF. Thus, 
in the FFF the normahzed, positive energy Dirac spinors of Eq. (57) take the rest- frame form 



4^\x) 








4^\x) 



1 





(60) 



It can be shown [13] that the general zeroth order solution to the Dirac equation in 
curved spacetime ipoi^) can be written as an amplitude times a normalized spinor ipo{x) — 
f{x)(po{x) (see also [22]), where f{x) satisfies the equation daf{x)u°'{x) = —9{x)f{x), 
and 9{x) = VQ'u"(a;) is the expansion of the cross section of a congruence ("tube") of 
timelike geodesies [25]. The normalized spinor fo{x) (i.e. ^q{x) (po{x) = 1 where <fo{x) = 
(f^{x)j^), satisfies u"{x)Da(po{x) — which states that (fo{x) is parallel propagated along the 
congruence. The results are most easily proved by introducing Riemann normal coordinates 
along an arbitrarily chosen timelike geodesic of the congruence along which the spin-| travels 
and for which the spinor connection ra{x) vanishes (see Eq.(40) and Eq.(41)). 

In this work, we will not introduce such local laboratory coordinates (which have limited 
spacetime ranges) and instead perform computations in the world coordinates x" of the 
metric ga/si^), and project world tensors to LIF tensors by means of the tetrad, as in Eq.(5). 
It is easy, and instructive to show that (po{x) is parallel transported along the congruence 
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in general coordinates. Taking ipo{x) — ipg (x) in the FFF from Eq.(60), with no loss in 
generality, we have 

u''ix)D^4^\x) = u''{x)[d^ + r^{x)]i;i^\x)^u''{x)r^{x)4^\x), 

= lu'^{x)uj^d^)a''4^\x) = '-Vace%{x)[u'^{x)Vaei-y'4'\^) = 0, 

where in the second and third equalities we have used the constancy of 4^\x) in the FFF, the 
definition of the spinor connection Taix) in Eq.(40) and the definition of the spin connection 
C(;^.g(x) in Eq.(21). The last equality follows from the definition of the FFF in Eq.(9), that 
the entire tetrad is parallel transported along the worldline of the FFF. Further, we note 
that from Eq.(23), x^i,{x) = in the FFF and additionally, w°(a;)r«(a;) = u^{x)ra{x) = 
e^°'{x)ra{x) — T^{x) — 0, while the spatial T'-{x) are in general non-zero, unless one uses 
locally adapted coordinates (e.g. Riemann normal coordinates). 

For the 0{h) solution to the DE, one notes that the operator in the square brackets in the 
non-homogeneous Eq.(50) is the same one that appears in the 0{1) homogeneous Eq.(49). In 
order for Eq.(50) to have a non-trivial solution ipi{x), all the solutions of the corresponding 
transposed homogeneous equation Eq. (49) must be orthogonal to the inhomogeneity on the 
right hand side of Eq.(50) (Fredholm alternative). This solvabihty condition for Eq.(50) 
becomes [13] 4^\x)j'^{x)DailJo{x) — 0, for a — {T, i}- In the following, we will only 
need the (9(1) solutions of Eq.(49) for calculating the 0{h) corrections the 4- velocity and 
4-acceleration of the spin-^ particle's trajectory. We therefore take as our solution to the 
DE in curved spacetime the wavefunction 

'iIj{x) = 'ijjl)'\x) 6xp 
with 4^\x) given by Eq.(60) in the FFF. 

B. Quantum corrections to the classical 4- velocity and 4-acceleration 

1. 4-velocity 

The designation of u°'{x) as the particle's 4- velocity arises from the equivalence principle 
for classical particles without spin. To exhibit the influence of the particle's quantum me- 
chanical spin on its trajectory one postulates that the motion of the spin-| particle through 
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- / p^,{x)dx^ 



a 



{T,U, i>{x)^l^{x) = 1. 



(61) 



the curved spacetime is determined by its conserved Dirac probability current j'^{x) [13] 

fix) = ip{x)r{xmx), D^rix) = o. (62) 

Performing a Gordon decomposition [26], the current can be written as the sum 

r{x)=j^{x)+jUx) (63) 

of a convection current i"(a;), and an internal magnetization current j^ix) defined by the 
minimal coupling curved spacetime generalization of their flat spacetime definitions 

Jcix) = ^ [^P{x)Do.^lJ{x) - (^^^(a;)) ^(x)] , (64) 

and 

JUx) = ^Hx)^"^{x)i:{x), a-^{x) = [r{x),^^{x)]. (65) 

We now let the convection current j^ix) define a congruence of timelike curves for the 
free motion of spin-| particles with tangent v^lx) 

j^{x) = v'^ix) - ^"(x) + 0{h), v''{x)va{x) = 1 + 0{h^). (66) 

The normalization of the 0{h) corrected v"{x) (right hand equality in Eq.(66)) will be re- 
quired for a self-consistent definition of a 4-velocity, and will be explicitly demonstrated 
in the examples considered in the subsequent sections. Using Eq.(61) and Daip{x) — 
D^i^^^^ {x) + {im/h) u^{x) 4""^ {x) ex.^{iS{x) /h) we obtain 

v^{x) = u^{x) + ^ [^^^\x)D^^l^^^\x) - {D^i^^^\x)) ^l^^^\x)] + 0(^1^). (67) 

We can simplify the above formula for Va{x) by noting that Dct'ip^\x) — {pa + 
V a{xy) 'il}^\x) — Va{x)'^^\x) and its adjoint Da,'4^^\x) — —'il)^\x)V ^{x) where we have 
used the constancy of iP^q\x) in the FFF, '4)^q\x) = i>l)'''\x)'y^ and 7°rJ^(a;)7° = —Taix), 
which yields 

va{x) = ua{x) + — ^a{x) 4''\x) + 0{h'). (68) 

mi 

Fig. (2) illustrates the spin-coupled non-geodesic trajectory of the particle with 4-velocity 
v(a;) perturbed away from its geodesic trajectory with 4-velocity u(x) when the spin of the 
particle is ignored. 
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FIG. 2: (Color online) The middle solid line represents the geodesic trajectory {a.{x) = 0) of 
the observer, riding along with the particle, with 4- velocity e^lx) = u(x) (black dashed arrow), 
if the spin of the particle is ignored (or considered as spin 0). The dashed red line represents 
the perturbed, non-geodesic trajectory (a(x) ^ 0) with tangent v(x) (solid red arrow), when 
the spin of the particle, which couples to the spacetime, is taken into account. The 4-velocity 
v(x) determines the particle's Dirac (convective) current ]c{x) = v(x) = e^^\x) which quantum 
mechanically defines the particle's motion through the spacetime. (Only the temporal axis of the 
C'(7i)-corrected FWF tetrad e^^\x), from which a null Wigner rotation is measured, is shown). 

We next compute the quantum mechanically corrected acceleration from the velocity 
correction Eq.(67) 

O-aix) = V^{x) DpVaix) = 2t>^ (x) f^] (x) , (69) 

which is the generalization of Eq.(8) when the 4-velocity contains both world vectors and 
spinors. The second equality in Eq.(69) follows from differentiating the normalization con- 
dition for f"(x) which allows one to add to the middle expression an identically zero term 
of the form v^{x)DaVf3{x) = 0. In differentiating Va{x), one must be mindful of the way the 
total covariant derivative Df3 acts on the individual vector Eq.(38) and spinor terms Eq.(39). 
A detailed derivation given in Appendix B reveals that the only nonzero terms that remain 
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upon differentiation of Va{x) are 

= -^Rap^six)uf{x)a^\x), (70) 

where we have made use of Eq.(44) for the commutator of total covariant derivative 
D[a, Dfsj = |[Da, D/}]. Equation (70) is the generahzed force equation for the v°(x) congru- 
ences and describes the deviation from the classical geodesic motion (described by u°'{x)) due 
to the coupling of the particle's spin to the curvature. The acceleration a°(x) = g°'^{x)afi[x) 
is force per unit mass that the spin-| particle experiences due to the coupling of its spin to 
its motion. Due to this non-zero, albeit small, acceleration the spin-| particle is no longer 
in a FFF (Eq.(9)) in which the tetrad is parallel transported along the particle's geodesic 
worldline (if one ignores the particle's spin), and is instead more appropriately described by 
a tetrad which is FW transported along the v"-{x) congruence by Eq.(lO). In the follow- 
ing sections we explore specific timelike worldlines in the Schwarzschild metric, which are 
geodesies if spin is ignored, and the quantum corrections to the velocity and acceleration 
given by Eq.(68) and Eq.(70). We explore the implications of the coupling of the particle's 
spin to the gravitational curvature on the Wigner rotation of the particle's spin and later 
its effect upon entanglement. 



VII. QUANTUM CORRECTIONS TO CIRCULAR GEODESIC MOTION IN THE 
SCHWARZSCHILD METRIC 

We now consider the quantum corrections to some specific classical geodesic motion (when 
spin is ignored) in the static spherically symmetric Schwarzschild metric 

ds" = (1 - 2M/r) de - ^ _ ^^^^ dr^ - {dO^ + sin^ dcj)^) (71) 

where x"' = {t,r,6,(f)), and = 2M = 2GM/c'^ is the Schwarzschild radius of the event 
horizon produced by the centrally located gravitating object of mass M (e.g. for the Earth 
Ts® = 0.89 cm, and for the Sun r^® = 2.96 km). Since the Schwarzschild metric is inde- 
pendent of the coordinates t (static) and (axial-symmetric) there exist two corresponding 
Killing vectors £ = (1, 0, 0, 0) and r] — (0, 0, 0, 1), respectively. These KiUing vectors can be 
used to define two integrals of the geodesic motion, by taking their scalar product with the 
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4- velocity u tangent to the geodesic trajectory 

e = ^ • M = (1 - 2M/r) 2M/r)-^, 

or 

I = -r^u^r'^ sin^ ^u"^ = sin^ 9 ^. (72) 

ar 

where u — (t**, li'', , u'^) in a coordinate basis. The quantity e in Eq.(72) is the total energy 
per unit rest mass, while / is the orbital angular momentum per unit rest mass. Both of 
these quantities remain constant along the geodesic trajectory. As in Newtonian central 
force motion, the constancy of the orbital angular momentum allows us to consider, without 
loss of generality, geodesic motion in the equatorial plane 9 — '!r/2 implying — dO/dr — 0. 
Writing and in terms of e and I respectively, and substituting into the equation for the 
normalization of the 4- velocity u ■ u — 1 produces a generalization of the Newtonian radial 
energy equation of the form (see Hartle in [16]) 



where the effective potential Veff{r) is given by 



V.,f(r) = I 



_M V_ _ Ml^ 



The first two terms of K//(^) in Eq.(74) are just the central gravitational potential and 
the centrifugal barrier potential of Newtonian mechanics. The third term in Eq.(74) rep- 
resents the general relativistic correction to Newtonian trajectories, which becomes more 
pronounced as r — > 2M. Extrema of Ve//(r) give rise to stable (+ sign) and unstable (- 
sign) circular orbits 

M / / V / ^ 12 
^circ = ^ 77 1 ± 777717? ■ (75) 



The innermost stable circular orbit (ISCO) occurs when (l/M)"^ = 12 or rjsco — 6M = Sr^. 



A. Circular geodesies and geodetic precession of gyroscopes 

For the case of circular geodesies r = R, using = dr/dr = in Eq.(73) along with 
Eq.(75) one can deduce the relations 

1-2M/R , VMR 
e = ' I = ^^^^=, on r = i?. 76) 

- 3M/R yi - 3M/R 
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We define the quantity Q = d(t)/dtas the orbital angular momentum of the orbit with respect 
to an observer at spatial infinity, where the proper time r (the observer's locally measured 
time) is equal to the Schwarzschild coordinate time t. We then have 

d(/> _ d(t>/dT _ l/R^ fW 

dt dt/dr e/{l-2M/R) M R^' ^ ^ 

Equation (77) is simply the statement of Kepler's third law that the square of the orbital 
period P — 27r/Q is proportional to the cube of the orbit's radius, oc r^. For a circular 

geodesic, the observer's 4-velocity takes the form u = (0, 0, m*^, m*), which from Eq.(76) and 
the normalization condition u ■ u = 1 yields 



I 00 ^ 

yi-3M/i?' ' ' y/1 - 3M/R J 



r=R 



where we have defined 



Note that for M = 0, the frequency Q'{r) contains the time dilation effect y 1 — V| = 
^/l — {rVl)^ for a particle executing circular motion with tangential velocity = rf2 (as 
measured from spatial infinity). 

For each circular orbit defined hy r — R simple algebra reveals that the covariant compo- 
nents of the 4-velocity are given by Ua — (e, 0, 0, 1) which by Eq.(55) shows that the action 
is given by S{x) — —{et-\-l(f)) appropriate for a circular orbit derived from the classical 
Hamilton- Jacobi equation Eq.(52). 

The interpretation of Vt'{r) can be inferred from solving for the spatial axes e^(a;) by 
treating the FFF condition Eq.(9) as a set of ordinary differential equations, with eQ(a;) 
given by Eq.(78). We consider the initial conditions such that at t = 0, er{x) points 
in the radial direction appropriate for a stationary observer at r = i?, = (for which 
^stat ^ _ 2M/i?)-i/2, 0, 0, 0)), i.e. ef^t = 0) = (0, (1 - 2M/Ry/'^, 0, 0). The solutions 
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to these equations yields the observer's tetrad 



e{x) 



ef{x) 











n'(r) 



Vl-2M/r n'{r) 




sin(Q'(r)0 ^1 - 2M/r cos{Q'{r) t) 

1 

r 



^ ^ sin(fi'(r) t) 



rf2'(r) 











cos{n'{r)t) ^yl - 2M/r sm{n'{r) t) 



y/1 - 2M/rVt 



rVt'ir) 



cos{n'{r) t) 



■sJl-2M/rQ.'{r) 

Prom the first row of Eq.(80) we observe that after one orbital period t — P — 27r/Q, the 
vector ef{P) does not return to its initial radial direction e^:(0). Rather, it is rotated in the 
direction of the orbital rotation by an angle 



'geod 



27r 



1 - 



271 



1 - yi - 3M/R 



(81) 



In fact, with respect to the initial, static orthonormal basis e{t — 0) (which is only a FFF 
frame when Jl'(r) t — 27rn) in which the spatial portions of the spatial basis vectors 6^(0), 
e^(0) and e^(0) point directly along the r, 9 and directions, respectively we have 

er{t) = cos(n'(r) t) er{0) - sin(n'(r) t) e^(0), 
e^{t) = sin(Q'(r) t) e^(0) + cos(Q'(r) t) e^(0). 

This effect is called (deSitter) geodetic precession. As discussed earlier, each of the spatial 
axes in Eq.(80) can be considered as a gyroscope defining a local axis in the observer's 
laboratory. Equation (81) implies that these axes precess by an amount A^^eod per orbit. 

It is important to note that we have written the tetrad components in Eq.(80) in terms 
of fl'{r) as a function of r given in Eq.(79), with Vl = ^JM/R^ treated constant. For a 
specific circular geodesic orbit, r and R are interchangable. However, when computing the 
quantum corrections to the velocity Eq.(67) and to the acceleration Eq.(69) we will need 
to compute the derivative of the tetrad, especially the variation in r. The expression for 
Q'(r) given in Eq.(79) has the important property that its r derivative evaluated on a given 
circular geodesic is zero, i.e. drO:'{r)\r=R — 0, upon using the definition of Q, given in Eq.(77). 
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The last point above can be more clearly understood by noting that the most general form 
for a stationary metric (i.e Qapix) indpendent of time, and gQi{x) 7^ in general, indicating 
a rotation of the spacetime) is (where we have temporarily reintroduced factors of c) 

ds'^ = e^*/"' (^c dt - dx^ - hj dx' dx^, i = (1, 2, 3), (82) 

in the rotating frame or lattice coordinates = (a;° = ct, X ^ X ^ X Here $ is analogous 
to the gravitational potential, but now includes length contraction due to the rotation. The 
spatial vector w — {wi,W2,W3) characterizes the rotation, and the 3x3 matrix kij is the 
spatial metric of the rotating spacetime. 

The gravitational field (acceleration) a experienced by an observer with fixed coordinates 
X = constant in the rotating frame has magnitude 

\a\ = [k'^{d,^){dj<^)Y^\ (83) 

where k^^ is the inverse matrix of the spatial metric kij. Further, if we put a gyroscope at 

— * 

X — constant in this rotating frame, it will precess at a proper angular velocity flgyro with 
magnitude 

relative to a freely falling frame (local inertial frame) [27]. The term in the square brackets 
is just the curved space generalization of the magnitude of the 3- vector curl w describing the 
rotation. In fact, for non-relativistic velocities we have the approximate expressions 

a — V$, ^gyro — TT curl'f/;, (85) 

and the total acceleration experienced by an observer at fixed lattice coordinates x is 

— * 

atotai ^ -a + 2v X Qgyro, (86) 

where — dx^ /dt. Thus, the total acceleration experience by an observer at fixed lattice 
coordinates x is the sum of inertial and Coriolis acceleration. In general, —a is the accel- 
eration and —^gyro is the rotation rate of the FFF with respect to the lattice (frame) with 
coordinates x". 

For the case of the Schwarzschild metric Eq.(71), we can define (restoring c = 1) 

= + d(j)^dip + Qdt, (87) 

31 



where Q is the constant angular velocity of a rotating lattice with coordinates x'°' — 
{t, r, 9, </?') with respect to the freely falling frame. After some straightforward algebra, 
the line element Eq.(71) takes the form (again for 9 — 71/2) 

ds^^il r^n^ ] \dt- z ^777 ^ V r-TTT ^7777 '-^ \ 



2M/r-r202 i-2M/r l-2M/r-r^n^ 

(88) 

From Eq.(83) the points of zero acceleration \a\, corresponding to dr^{r) = correspond 
to a geodesic, i.e. a free particle can remain at rest there (in the new coordinates) and the 
worldline of the lattice point is a circular geodesic. A simple calculation yields 

dMr)\r=R^O ^ = (89) 
A calculation of flgyro = l^ayrol from Eq.(84) yields 

^gyro (90) 

The above is only an apparent coincidence since ^^gyro is a proper rotation rate relative 
to the FFF, while is a coordinate rotation rate (with respect to t). At fixed lattice 
coordinate .-f' on the circular geodesic r = Rwe have At = (1 — SM/RY^"^ At = dt/'f where 
7~^(r) = (1 — 2M/r — r^O^). Therefore, after one orbital period At = 2tt/Q,, a gyroscope 
at fixed lattice point x' traverses an angle 

cpgeo - O Ar = — = n^/l - 3M/R ^ = 27rVl - 3M/R. (91) 

The precession of the gyroscope per orbital revolution relative to original lattice is then 
^fgeo — 2% — (figeo which is the same as Eq.(81). 



B. Velocity and acceleration corrections 

To compute the spin-orbit quantum correction to the covariant components of the velocity 

Va{x) = Ua{x) +h5Va{x), SVa{x) = 'lpo{x) Taix) 1po{x) + 0{h^), (92) 

im 

we use Eq.(40) for the spinor connection Taix), which in turn uses the spin connection 
a; -^(x) computed from the definition Eq.(40) utilizing Eq.(34) and Eq.(35). For 'ijjQ^\x) 
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we form the most general constant component spinor from the positive energy solutions for 
spin up and spin down Eq.(60), in the particle's FFF 

/ cosC/2 \ 



\ 



sinC/2 





(93) 



/ 



where (C, ^) arc the constant polar and azimuthal angles relative to the quantization axis. In 
the observer's rest frame the spin points along the direction n — (sin ^ cos </?, sin C, sin </?, cos C,) . 
For motion in the equatorial plane 9 — Tr/2,we choose the quantization axis, the local z-axis, 
to be perpendicular to the orbital plane along — e^, and hence the local x and y axes to be 
along ef and respectively. 

With the spinor Eq.(93), the quantum velocity correction 6va{x) has the general form 



5Va{x) 



1 

i m 



cos2(C/2) rl^x) + sin^(C/2) ri^{x) + - sin C {e'^ FI^^) + e''^ r^J{x)) 



(94) 

where V^^lx) = ip^ \x) Ta{x) '4)^\x), with il)\^\x) for a', c G the spin up and spin 

down spinors from Eq.(60). For the specific case of the Schwarzschild metric Eq.(71), the 
spinor connection has the properties 



Ti\x) 



0, 



Fit, 



X] 



0, 



r^^ [x) 



Scwharzschild metric 



Ti\x) = FjT(x) = 0, T'^{x)^-[Tl\x)y . 



A lengthy, but straightforward calculation yields the correction 5va{x) for arbitrary radius 



5va{x) 



5Vr{x) 

5ve{x) 



( 



QR cos C 



2m^yl - 3M/R 



- 2M/R sm{Q't - (p) sin C 
2mi?2 



(1 -2M/i?) cosC 



(95) 
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where Q = -yjM/R^ and = Q (1 - SM/RY^'^. 

In general, the velocity correction in Eq.(95) depends on the orientation of the spin in the 
FFF. For the case of spin up {( = 0,(p = 0) and spin down {( = TT,(p = 0), 5va{x) has only 
SV(f,{x) and Svt{x) components, and thus the motion is along the circular orbit. For values 
of < C < TT, there is a non-zero ()vg{x) component, corresponding to small oscillation in 
the local z direction at frequency fl'{R). 

A straightforward calculation of the acceleration aa{x) proceeds from Eq.(70). Here the 
Riemann curvature tensor Rapysix) is computed from the Christoffel connection directly 
from the metric from Eq.(46) and Eq.(7), and the world Dirac matrices a^^{x) are calculated 
from the usual Dirac matrices in fiat Minkowski spacetime a'^'^ from Eq.(65) and Eq.(42). 
This yields 



aa{x) 



I (^t{x) \ 
ar{x) 
ae{x) 



( 



\ 



SQ^R cos C 



V 



2m(l-3M/R) 

Vt^/l^^2M/R cos(Q' t - ip) sin C 
2m 





(96) 



/ 



For pure spin up or spin down in the FFF, the acceleration is strictly in the radial direction, 
corresponding to the purely circular velocity corrections discussed above. The acceleration 
aa{x) is the force per unit mass that the spin-| particle experiences due to the coupling 
of its spin to the gravitational curvature. This coupling produces a change of the particle 
motion from the geodesic to to non- geodesic to 0{h). With respect to the FFF the 
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acceleration has the components (a = a"(x) ea(x)) on r = i? 

/ a*(x) \ 



a (x) 



\ 



n 



J 



Sn^R y/l - 2M/R cos{n' t) COS ( 
2m{l-3M/R) 



- 2M/R cos(l^' t - ip) sin C 



2mR 



Sn^R y/1 - 2M/R sin(Q' t) cos C 



V 



(97) 



2m{l-3M/R) 

As a useful consistency check, it is worth noting that aa{x) arises directly from the 
differentiation of h6va{x) in Eq.(92) keeping in mind the action of the covariant derivative 
on the different type tensorial/spinor quantities. Thus from Eq.(92) we have 



h 



hDf,6v^{x) = — Df,{iji;'\x) r„(x) <^(a;) 
mt ^ 

h 



= — ( {Dpi^^^\x)) V^{x) 4''\x) + VpT^ix) 4^\x) + 4^\x) r^{x) D^4^\x)) , 

( - Tp{x)) T^{x) 4''\x) + #^(0;) V^r„(a;) 4''\x) + #^(0;) T^ix) Tfs{x)4-\ 



rrn 

h_ 

mi 

h 



= ^#(^) i^p^aix) + [r^ix), r^ix)]) V'r(x), 

so that upon anti-symmetrization on the indices a and (3 we have 

a„(a;) = u''{x)2D^6v^] = 1^ {x) 4''\x) (V[^r,](a;) + [T^{x),T^{x)]) 4''\x). (99) 

mi 

In the above Dp acting on the spinor connection Taix) is just the Riemann covariant deriva- 
tive as in Eq.(31), while its action on the spinor 4'^\x) is given by Eq.(39). The term in 
the parenthesis is just the explicit expression for 

[Dp,D^]4''\x) = '-Rp^,s{r{x))a^\x)4''\x), 

given in Eq.(B6) and hence yields the expression for aa{x) in Eq.(70) (note Rpa-yS — —Ra/Sys)- 



C. Corrections to the tetrad 

To find the 0(h) corrected tetrad ef'\x) that defines the FWF (instantaneous, non- 
rotating rest frame of the particle) in which the observer detects a null Wigner rotation. 
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one solves the FW transport equations Eq.(lO), setting eQ(x) = v{x) given by Eq.(92) and 
Eq.(93). Let s be any one of the three spatial tetrad vectors ej(x) with the property s-v — 0. 
We employ a Lindested-Poincare perturbation approach [28] which allows for a nonlinear 
frequency correction. The 4-vector s is expanded as s(x) = Sq{x) + hsi{x) + 0{h'^), and 
we define a new renormalized time ^ = ujt = (1 + nui + 0{Ti^)) T, with r the proper time. 
The 0{1) equations for Sq{x) reproduce the FFF tetrad of Eq.(80) where the argument of 
the sinusoidal functions are Q^, reducing to Qr in the limit cui ^ (which equals D,'{r)t 
from the discussion before Eq.(91)). 

The condition s ■ v — allows one to write Sq{x) in terms of Sq{x) and s\{x) in terms 
of both sf{x) and Sq{x) (for components written with respect to a coordinate basis). For 
motion in the equatorial plane 9 = 11/2 (so that Sq{x) = s^{x) = 0) the remaining equations 
for s'[{x) and (x) can be used to construct a second order equation for sl{x) which takes 
the form 



dlslix) 



+ n s\{x) = f{v,u;i)sl{x), 



where f{v,u!i) is independent of the time ^, and evaluated on the circular orbit r = R. 
The trivial particular solution Si{x) = can be found by using cui to make f{v,u!i) — 0. 
A detailed calculation reveals that for an arbitrary circular orbit of radius r and with v 
computed for a pure spin up spinor Eq.(60), 

^iW = 7777^ l + 7^h (100) 



4M Q'(r) V ^'(r) 

which, in addition, has the property that drUJi{r)\r=R = 0. Thus, the first order uniform 
expansion solution (see [28], pl25-126) for the FWF tetrad is the FFF tetrad of Eq.(80) 
with the frequency fl'{r) in the sinusoidal functions replaced by the renormalized frequency 

$'(r) = (1 + ;ia;i(r)) 0'(r). (101) 

The above implies that the geodetic precession rate per orbit Eq.(91) is increased to 



= 27r (1 + hu.ir)) Vl - 3M/i? (102) 

on a circular orbit of radius r = R. 

The calculation of the spatial tetrad axes e?'^ (x) , proceeds similarly, but yields comph- 
cated expressions. Most importantly, the correction to the geodetic precession of the spatial 
axes depends on the particular orientation n(C, ip) of the spin in the observer's local frame. 
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As stated at the beginning of this section, the observer would detect a zero Wigner rotation 
in the non-zero acceleration FWF ef'\x). 



D. Wigner rotation 

As discussed in the previous section and at the end of Section IV., an observer in the 
FWF will detect no Wigner rotation of the particle's spin. It is instructive to compute the 
the Wigner rotation as observed from the FFF in which the classical general relativistic 
motion is a geodesic when the particle's spin is ignored. The non-trivial portion of the 
Wigner rotation matrix '&\{x) Eq.(26) is calculated using p°'{x) — mv"'{x) — me"' ^{x) v^{x) 
and a"{x) = me"- ^{x) a^^{x) where v^^{x) = u^(x) + h5v^{x) and a^{x) = hSa^{x) are the 
0{h) corrected 4-velocity and 4-accclcration from Eq.(92), Eq.(93) and Eq.(94), and e"^(a:) 
is the transposed inverse of the FFF tetrad e.'^(a;) Eq.(80), of the observer. The expression 
for x'^i^i^) non-trivial portion of the infinitesimal LT A"-(a;) in Eq.(24) is obtained 

from Eq.(23) by replacing u{x) by v{x), i.e. x"^(a;) = — e"^(x) Vv e- '^(x). 

Keeping terms to 0(h), a straightforward calculation of Eq.(24) yields 



A«^(x) = [a'(x)vi(x)-v\x)ai,(x)]+x'i(x), 



= h 



{Sa%x) ui^ix) - u%x) 5ai^{x)) - e"^(a;) 5v^{x) ^(x) 



hSX".{x), 



(103) 



where we have made use the FFF condition Eq.(9) for the tetrad. The relevant portion of 
the infinitesimal Wigner transformation from Eq.(26) becomes 



X '.ix) H — 



= h 



v^ix) + 1 

xxi r ^ , SX'^{x)u.{x)-u'{x)SX.f^{x) 
dX '.{x) H — 



(104) 



u^{x) + 1 

In examining the terms in Eq.(102), we note that since the observer is using the FFF in 
which eQ{x) — u{x), the local components are simply u"'{x) — 5"^. By the discussion at the 
end of Section IV all the terms involving u^{x) in A*-.(x) and (x) vanish and we are left 
with 



X\{x)=x){x)=i}){x), 
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with x'-.(x) given by the second term in Eq.(lOl) involving 5v^{x) yielding 



( cos(Q't) sin(Q't) sinC -cosC ^ 

- cos(Q' t) sin(Q' t) sin C - %\^{^ t) sin C 

^ cosC sin^(r2't) sinC j 

, , (1-2M/R) , , 

In Eq.(105), x^(x) is the single, non-zero value X'i{^^ — ~X^3(^)) (where {1,2,3} = 
{x, — z,y} <-> {er,eg,e^, }) computed for a pure spin up spinor (C = 0), Eq.(60). ^{x) 
represents a spatial rotation about the local z-axis (— e^-axis) perpendicular to the plane of 
the orbit, of the FFF with respect to the FWF. 

For the orientation we have chosen above for the observer's local axes in the equatorial 
plane the infinitesimal spinor rotation matrix Eq.(27) becomes 



where 




\=^^^ ^^2 = - I ^ ^ j = -<^z, = 1^ . ^ \=<^y- (107) 
Here we have defined 

{ cosC/2 \ , . 

V'o= . \ ^ \v\x),a) (108) 

y e*'^ sinC/2 j 

as the upper two components of ■0o of Eq.(93) and which we can associate with the positive 
energy state |p'(x), a). In Eq.(106) we have indicated that D^^J^^ depends upon the Wigner 
rotation matrix W{x), which in turns depends on the zeroth-order spinor wave function ■j/'o 
(which is parallel transported along the circular geodesic with 4- velocity u[xy). For example, 
if we consider a pure spin up or pure spin down state, (C = 0, = 0) or (C = tt, = 0) 
respectively in Eq.(108), only 'Q^i — x^(i?) cos^ is non-zero. In these particular cases, 
is independent of time and hence the time ordering in Eq.(29) can be ignored, and the spinor 
rotation matrix can be integrated for finite rotations 

I, (109) 
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where the upper and lower signs are associated with spin up and spin down respectively, 
and r(i?, t) — Q'{R) t. However, for a general spin orientation < C < tt in the FFF Eq.(93) 
(e.g. (C = 7r/2,(y9 = 0) corresponds to the spinor pointing along e^, {( = 7r/2,(y9 = 7r/2) 
corresponds to the spinor pointing along e^), 'dij{x) is time dependent Eq.(105), and 
D^y^^ (W{x)] ipo) is not simply a rotation about the local z-axis, and can only be computed in- 
finitesimally. The state \p\x'),a) = U{A{x)) \p\x),a) = D%^\w{x'y,ipo) \p\x'),a') 
represents the Wigner rotation of the state \p^{x),a) as measured from the FFF at the point 
x' due to the 0{h) correction of its spin coupling to its motion. 

E. Entangled States 

Terashima and Ueda [11] considered the Wigner rotation of a spin-singlet state (with the 
local quantization axis along e^) created at given point on a non-geodesic equatorial circle, as 
one of the particles of the bipartite state circulated the orbit clockwise, and the other counter- 
clockwise. In this case the observers are taken to be a set of stationary observers situated 
around the orbit at r = i? (the only non-zero tetrad components are e^^(x) = (1 — 2M/r)^/^, 
e^^{x) = 1/r, ^^{x) = l/r, and e^^{x) = (1 — 2M/r)~^/^). Since the circular orbit was non- 
geodesic, a represented the external, non-gravitational 4-acceleration required to keep the 
particles on the orbit, with each particle experiencing equal and opposite spatial tangential 
velocities. Given the velocity and acceleration of each particle, they calculate the appropriate 
Wigner rotation and corresponding unitary transformation for each particle and then apply 
Di{W{^)) (8) D2{W{-^)) to the bipartite state. Here the subscripts 1 and 2 refer to the 
first and second particle comprising the bipartite state, which traverse the orbit from the 
point of origin = to = $(i?) and (f) = — $(i?) respectively. 

The result of this calculation is that the spin-singlet state is mixed with the spin-triplet 
state, and hence spin measurements in the same direction are not always anti-correlated in 
the local inertial frames at = ±$ (i.e. along the local axis e^). The general relativistic 
effects deteriorates the perfect anti-correlation in the directions that would be the same as 
each other if the spacetime were flat [11]. This deterioration is a consequence of the fact 
that in their scenario 

A^.(a;)7^X^(^)7^^y^), 
which results from the imposed accelerations and particular choice of stationary observers 
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(tetrads) situated around the orbit. 

For the problem that we consider, the couphng of the particle's spin to its motion produces 
a deviation of its motion from geodesic to 0{1), to non-geodesic to 0{h). Let us consider 
the case of an entangled bipartite state in which to 0{1) the individual particles traverse 
circular geodesies at slightly different radii r = R± 6R in the same direction, and then ask 
how the 0{h) corrections to the particles' orbit effect the bipartite state as observed from a 
co-circulating FFF situated at r = i? (see Fig. (3)). For our entangled state we consider the 




FIG. 3: (Color online) Bipartite entangled state |^'(x4., j;_)) located on equatorial circular geodesies 
r = Rziz AR (coordinates x± in text) when the spin of the particle is ignored. We consider the 
0{fi) spin-orbit coupling of the particle leading to non-geodesic motion on x± (indicated by dotted 
red lines), and we consider the Wigner rotation from the geodesic FFF at r = R. The dashed 
sinusoidal blue curve, perpendicular to the equatorial plane, is the non-geodesic motion when the 
spins of the particles are not along the local z-axis (e^). When the spins are oriented parallel or 
anti-parallel to the z-axis, the 0{h) velocity corrections are parallel or anti-parallel to the circular 
orbit. 
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following normalized bipartite state 

= cos(e/2) |p(x+),t) + e^* sin(e/2) \p{x+),i) |p(^-),T)- (HO) 

In Eq.(llO), p{x) represents the spatial portion p'^{x) of the local momentum j9"(x) with 
respect to the local tetrad e(a;). The argument x± indicates the circular orbits of radii 
r = R±SR, respectively. We have also included the arbitrary, constant angles (O, $) in the 
definition of our bipartite state. For (© = 0, $ = 7r/2) Eq.(llO) is a pure spin-singlet state 
(J = 0,mj = 0), and (© = 0, $ = 0) represents the triplet state (J = l,mj = 0), while 
for < < TT it is a linear combination of the spin-singlet and spin-triplet states of zero 
magnetic quantum number. 

For each particle we define the quantum correction to the velocity as a generalization of 
Eq.(64) for a multi-particle state 

n 



V, (X = 



: [^{x+,x_) D^{x+) (g) I-^{x+,x_) - {Da{x+)(^I'^{x+,x_)) '^{x+,x_)] , 



' 2mi 

V2{x) = ^ [^{X+,X^) I ® Da{X-)^{x+,X-) - {I ® Da{x^)^{x+,X^))^{X+,X^)] 



(111) 



This leads to expressions analogous to single particle velocity corrections Eq.(68) 

h 



Vla(x+) = Ua(x+) H : -^{x+,X-) (Ta(x+) l) '^(x+,X^), 

rrn \ / 

V2a{X-) = Ua{X-) + — ^{x+, X-) (l ® Ta{X-)) '^{x+, X-), (112) 



where it is important to note that the entire entangled wave function ^(a;+, x^) Eq.(llO) is 
used to calculate these local velocity corrections. An explicit calculation yields 

vi^{x+) = u^{x+) + — [cos^(e/2)rIT(x+) + sin2(e/2)rii(x+)] =«,(x+) + — coserIT(x+), 

rrn mi 

V2a{x-) = u^{x^) + — [cos2(e/2)riH2;-) + sin2(e/2)rTT(x_)] =K,(2;-)- — coserlT(x_), 

mi mi 

(113) 

where we have used T\^{x±) — —T]^{x±) for the Schwarzschild metric (section VII. B). 
Relative to the circular orbit at r = (which we will indicate in subsequent expressions 
with the argument {x) with no subscripts), the 0{K) velocity corrections 

5vJx±) = ±— cos e V^^{x±) (114) 
mi 
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are parallel or anti-parallel to the direction of motion for the two particles, due to the anti- 
correlation of the spins in the bipartite state Eq.(llO). Note that for a pure spin-singlet 
or spin-triplet entangled state (0 — 7r/2, with $ = 0, $ = 7r/2 respectively) the velocity 

correction is zero to 0{h), which implies no 0{h) Wigner rotation correction. 

Relative to the FFF for each circular orbit with tangent Ua{x±), the local 4- velocity 
components take the form u°'{x±) — S"'^. Thus, the results of section VII. D apply to each 
orbit x± and hence 

From Eq.(103) x*-(a;±) takes the reduced form 



(115) 



Since rlj{x±) is an expectation computed with a pure spin up state, again the only non- 
zero value of the x*-(^±) matrix is x\{^±)\c=o — —cosQx^{^±) Eq.(105), now evaluated at 



r — R± 5R. Let us define 



■d{x±) = '(?3i(x±) = X3i(^±) = cosGxU^i)- 



(116) 



As a result of the alternating signs in Sva{x±) Eq.(114), the Wigner rotation angle 'i?(a;±) is 
opposite for the two orbits x±. 

We now wish to describe the Wigner rotation of the bipartite state as seen from an 
observer in the FFF at r = with 4- velocity u°'{x). Taking 5R/R <^ 1 we can expand the 
terms in Eq.(114) and Eq.(115) to 0{5R) as 



dr 



to obtain 



X^Xi) = ±x){x) + SRAxY^), 

x'-Ax) = — ^ cose 

^ mi 

AxVx) = -—cose 
■' mi 



{x) V^^\x) Vpe.^{x) 



9 i , . 



T^^\x)V,e.^{x) + e\{x) [^^T^^\x)]V,e.^{x) 



d 



+ e\{x)T^^\x)^^ (V,e/(a;)) 
where Vl^^^x) = g^'^{x) T]^{x). Thus, Eq.(116) becomes 



(117) 



i?(x±) = ±^{x)^M{x), ^{x) = cosex^(^), S^{x) = SR Ax3i{x) (118) 
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In Eq.(117) and Eq.(118) the leading ± in front of X^ji^) "^i^) arises directly from the 
alternating sign in Eq.(114) for the particles on orbits of radii r — R± 6R. The additional 
± sign, resulting from expansion of the argument x±, combines with former ± to produce 
the 0{6R) with the same sign for both particles. 

Prom the FFF at radius r — R the unitary transformation of the bipartite state Eq.(llO) 

is 

D[W{'d{x) + A^{x))] ® D[W{ - -dix) + Ai?(a;))] = 

g gi/2(i?(x)+Ai?(a;))(iT 

The action of Eq.(119) upon the bipartite state Eq.(120) produces the Wigner rotated state 

\^'{x)) = cos(e/2) e-'^(^)'^" \p{x), t) \p{x), D + e"" sin(e/2) e*''^^)'^^ \p{x), i) \p{x), t), (120) 

as observed from a FFF observer circulating the geodesic circular orbit at r = i? with 4- 
velocity u°'{x). Note that to lowest order the 0{6R) have cancelled identically. Further, for 
the pure spin-singlet (6 = 7r/2,$ = 0) and pure spin-triplet state (6 = 7r/2,<l> = 7r/2) of 
zero magnetic quantum number, there is no observed Wigner rotation, since from Eq.(118) 

= cosex^(a;) |(e=7r/2) = 0. 
For any other observer with tetrad e'(a;), instantaneously coincident at r = i? with the 
FFF observer with tetrad e(,T), the observed Wigner rotation of the bipartite state is much 
more complicated. It can be obtained by finding the components of the 4-vclocity v{x) 
and 4-acceleration a{x) for each particle relative to the new observer's tetrad, which is 
calculated by a local (position dependent) LT between the two instantaneously coincident 
frames, e^(x) = A-''(x) e^(x). Finally, by using the first lines of Eq.(103) and Eq.(104) one 
can work out the Wigner rotation matrix as observed from the new frame. 



^-i/2[-i}(x)+A&{x))dT Q 

Q gi/2(— !9(x)+Ai?(x))dT 



• (119) 
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F. Comments on Spin-Momentum Entanglement and Wigner Rotation 



A general two particle state takes the form 



dp 




^ J J ^^^^^^ 

iPuP2)\pl{Xl),ai) \P2{X2),(T2), 
^ dpidp2 \9na2{Pl,P2)\'^ = 1, 



(121) 



where gaia2{Pi-iP2) is the joint spin-momentum distribution function and dp is the (local) 
Lorentz invariant integration measure. Prom Eq.(14), an infinitesimal LLT U — U{K{xi)) ® 
U{K{x2)) of the state \^{xi,X2)) will mix spin a and 4-momentum p since the Wigner 
rotation angle is momentum dependent, the same as the flat spacetime result discussed by 
Gingrich and Adami [6]. As a result, the reduced two particle spin density matrix formed by 
tracing out the momentum of the particles will exhibit an observer (Lorentz transformation) 
dependent entanglement, measured e.g. by Wootter's concurrence. 

The results of the previous sections indicate that, due to the spin-curvature coupling of 
the massive particle in CST resulting from the particle's orbit (4- velocity) being determined 
from its Dirac current, the Wigner rotation also depends on the initial orientation of the 
particle's spin in its local frame (sec the wavefunction Eq.(93), the velocity Eq.(95) and 
acceleration Eq.(96) corrections, and subsequent expressions). Thus, in CST the observed 
Wigner rotation is a function of not only of the LLT A and 4-momentum p, but also the 
spin orientation fi(C, </?) in the local (laboratory) frame, which we indicate as 



VIII. RADIALLY INFALLING GEODESIC MOTION IN THE 

SCHWARZSCHILD METRIC 

To explicitly demonstrate the local spin orientation dependence of the Wigner rotation 
Eq.(122), we consider the simpler case of the quantum corrections to the FFF of a radially 
infalling particle, and the 0{h) corrected tetrad in which the Wigner rotation would be mea- 
sured null. For a particle dropped in from spatial infinity in the equatorial plane {9 — 7r/2) 



W = W{A,p,n), 



in CST. 



(122) 
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with zero velocity and zero orbital angular momentum {e — 1,1 — 0), the radial energy equa- 
tion Eq.(73) with the ansatz u = {u\ u'^, 0, 0) yields u = (1/(1 - 2M/r), -(2M/r)V2, Q, 0), 
where the minus sign in indicates a radially infalling geodesic. The freely falling frame 
tetrad (FFF) in the equatorial plane is given by 

1 



e{x) 



2M/r 



^2M/r 
' 1 - 2M/r 







-^/2Mji 











1/r 



e^{x) 



ei^x) 



e^[x) 



63(2;) 



(123) 



1/r 

which satisfies the FFF condition Eq.(9), VvL^a{x) = 0, where ^^{x) = u{x), and the 
orthonormalization condition e{x) ■ g{x) ■ e^{x) = ij. (Note, we have used the same local 
axes designation as in Fig. (3)). 

Using the general spin orientation in the FFF given by Eq.(93), corresponding to a 
local spin direction ofn — (sin C cos </?, sin C sin (p, cos Q , we obtain velocity correction Sv'^lx) 
Eq.(94) and the acceleration correction 6a°'{x) from Eq.(99) 

/ vt{x) \ ( 1/(1 -2M/r) \ 



v'^{x) = u'^{x)^h5v'^{x) 



Vr{x) 

ve{x) 



-V2M/r 

~h sin C sin 99/ (4Mr^) 
\^ -^cosC/(4Mr^) ) 



8a'^{x) = 0. 



(124) 



For a spin orientation in the FFF in the local x-z {f-6) plane in which (/? = 0, the effect of 
the velocity correction is a deflection along the local y (0) direction in the equatorial plane 
by 5v^{x) — — cosC/(4Mr^). For sin C sine/? :^{){ny^ 0) the velocity correction takes the 
motion of the particle out of the equatorial plane. For pure spin up/spin down {C, — (0, tt), 
respectively) we have 5vY{x) — ^h/{4:Mr'^), as depicted in Fig.(4). Since the acceleration 
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FIG. 4: (Color online) Bipartite entangled singlet state |^'(a;,x))i2 Eq.(127) on a radially infalling 
geodesic with FFF e(x) with 4- velocity e^^x) = u{x) Eq.(123). The coupling of the spin to the 
spacetime curvature creates a spin dependent deflection of the particles with the new 4-velocity 
(Dirac current) v{x) = u[x) + fi6v[x), with 6v^{x) = — cosC/(4Mr^) Eq.(124), which is along the 
local =Fy-axis for pure spin up (C = 0, (/9 = 0) and spin down {( = Tr,ip = 0) orientations (where 
n = (sin C cos 99, sin ^ sine/?, cos C) is the spin orientation in the local inertial frame at x). Under a 
LLT, the state ?7(A)|^'(x, x))i2 = \^{x'_^_, x'_))i2 at the transformed points x x'^. = x + v dr 
remains a spin singlet state if the two observers use different C'(?i)-corrected FFF tetrads (since 
5a{x) = 0), e(^)(x'_; C = 0, 99 = 0), for spin up at x-, and e!^^\x'j^] (" = vr, 99 = 0) for spin down at 

Xj^. 

correction is zero in Eq.(124), the Wigner rotation is of the form given in Eq.(115) for an 
arbitrary tetrad e(x), and arises solely from the rotation of the tetrad x^j{^) ^i^nd the velocity 
correction 5v"'{x). Since 5a°(x) = 0, we can solve for an 0{fi) corrected tetrad e^^)(a;) which 
is the FFF in which the observer measures a zero Wigner rotation. Solving 

FFF in which W = : Vv &f\x) = 0, ef\x) = v{x), (125) 
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yields 



ef^(x) 



ef\x) 



1 - 2M/r 

_ V2M/r 
~ 1 - 2M/r 

sin ( sin (/? 
4Mr(l - 2M/r) 

cosC 

L"4Mr(l-2M/r) 



. sinC sin 09 , cos^ 
-n — „ — —a 



^ sin C sin cp ^ 



h 



-h (r - M) 



sin( simp 



1/r 



4Mr2 
cosC 



n 



cosC 



which satisfies the normahzation condition e^'^^ I 



X ■ K Uc ■ e 



{n)T, 



1/r 
(126) 

x) = r] + 0{h^). In Eq.(125) 



we see expUcitly that the FFF e^^\x') = e^^\x';C,ip) in which W = at the transformed 
point x' depends on the spin orientation n{(, ip) at x. 

An imphcation for entangled states is as follows. Consider, for example, the spin singlet 
state at the spacetime point x (see Eq.(llO)) 



|*(X, X))i2 = \p{x), t)l \p{x), i)2 - \p{x), i)i \p{x), t)2- 



(127) 



Since 6vY{x) = ^l/(4Mr2), the spin up states will be deflected along the local — y-axis, 
while the spin up particles will be deflected along the local +y-axis. In order for the state 
to be observed with zero Wigner rotation, i.e. as the state 



|*«,x'_))i2 = bV-),t)i |P«)4)2 - bV+)4)i bV-),t)2. 



(128) 



where x'^ indicates the transformed particle positions with small displacements along the 
local ±y-axis, we need the two local Lorentz observers (or local inertial frames, LIF) to be 
using different tetrads e*^^^ at x'_^, i.e. the LIF at x^ using the tetrad appropriate for spin 
up e^^\x'_-X — 0,(p — 0), and the LIF at x+ using the tetrad appropriate for spin down 
e'^^\x'_^_-X — n,(p — 0), see Fig. (4). We cannot use a single class of observers throughout 
spacetime, since e^^\x') depends on the initial spin orientation n{(,(p) at x. For a pair 
of observers other than e^^^\x';(^ = {vr, 0},(/9 = 0) at x'^, the general LLT Eq.(28a) and 
Eq.(28b) indicates that transformed state U{A{x)) \'i/{x,x))i2 would be an entangled state 
composed primarily of a singlet state, with an 0{h) superposition of a triplet state. Such 
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observer's would detect an apparent decrease in the maximum EPR correlation for the two 
particles at while those using e^^\x'; ( — {tt, 0}, (p — 0) would not. 

IX. SUMMARY AND CONCLUSIONS 

In flat spacetime the positive energy, single particle state \p, a) of a massive particle is 
given by its momentum, and spin-j components a along some quantization axis. Under a LT 
A, describing the same flat spacetime from an inertial reference frame moving with constant 
velocity (zero acceleration) relative to the original inertial frame, the state is transformed 
by Eq.(ll), in which W{A,p) is the Wigner rotation angle of the {2j + 1) x {2j + 1) 0(3) 
rotation matrix, which mixes up the spin components of the transformed state. If we consider 
wavepackct states composed of an integration of single particle states over a momentum 
distribution, then the momentum dependent Wigner rotation angle W{A,p) will transform 
each momentum component differently. For a general wavepacket comprised of a joint 
distribution over two particle states, the unitary transformation in Eq.(ll) will lead to a 
spin-momentum entanglement, such that if we consider the reduced spin density matrix of 
a bipartite state, there will be degradation of spin-spin entanglement. 

In curved spacetime, we replace the positive energy, single particle state of fiat spacetime 
with the local state \'p{x),a) valid in the locally fiat Lorentz tangent plane to the the CST 
at the point x. By the equivalence principle, the laws of SR hold in this tangent plane. 
Measurements of properties of this state are made from the reference frame of a massive 
observer with 4- velocity Ugbsix), instantaneously collocated at the spacetime point x. We 
define this observer by the orthonormal tetrad ea{x) which comprise the four axes of his 
local laboratory. The three axes G:-{x) comprise the spatial axes at the origin of the ob- 
server's local laboratory, and the temporal axis GqIx) = Uobs{x), tangent to the observer's 
worldline, determines the local rate at which his clock ticks. The observer measures the 
local components p"(x) of the 4- momentum p{x) — mu(x) of a particle passing through his 
local laboratory at x by projecting p{x) onto the tetrad via p"(x) = ^{x) p°'{x), where 
E{x) — p^{x) is energy of the passing particle as measured by the observer, and p^{x) are 
locally measured 3-momcntum components. Under a local Lorentz transformation A(a;), 
which transforms between observers in different states of motion collocated at x (i.e. sta- 
tionary, freely falling, circular orbit, or under arbitrary acceleration), the state \p{x),a) 
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transforms under a local Wigner rotation W{A{x),p{x)) Eq.(14), which generahzes the flat 
spacetime result of Eq.(ll) by the equivalence principle. The same SR considerations of 
spin-momentum entanglement for wavepackets states hold as well in CST, but now locally 
at each spacetime point x. 

In this work we have considered massive particles of spin i. Since the observed particle 
is massive, we can always find a non-rotating, (in general) accelerating observer that is 
instantaneously at rest with the particle, called a Fermi- Walker frame (FWF). If the particle 
follows a geodesic (zero acceleration) the FWF reduces to the particle's freely falling frame 
(FFF). We have shown that in the FWF, the Wigner rotation angle is null. In any other 
frame, the observer would measure a non-zero Wigner rotation for the spin of the particle. 

It is a postulate of General Relativity that the force free motion of particles are geodesies, 
i.e. trajectories for which the acceleration is zero. This is true if the particle is assumed to 
possess no spin. If the particle does possess spin (even classically), the spin of the particle 
couples to the curvature and creates non-geodesic motion. In this work, we have considered 
the motion of a quantum spin | particle, as determined by its Dirac current, leading to non- 
geodesic motion to first order in h. We have found that the momentum (and acceleration) 
corrections depend on the initial spin orientation of the particle in its local frame, and hence 
the momentum dependent Wigner rotation will transform each of these spin orientations 
differently. Whereas in flat spacetime the Wigner rotation angle W depends on both the 
Lorentz transformation A and the particle's momentum p, in CST W depends on the local 
analogues of these quantities, as well as on the spin orientation n of the particle in its local 
frame. In addition, we have explored the evolution of the entanglement of a bipartite state 
on two infinitesimally close circular orbits, as observed from the FFF of circular orbit at the 
averaged radius. We have also shown that an entangled Bell state formed at the spacetime 
point X for a pair of collocated freely falling observers in the equatorial Schwarzschild plane, 
will be observed to be in the same entangled state, i.e. observed to have zero Wigner rotation 
under a local Lorentz transformation, if at the spatially separated transformed spacetime 
points, it is observed by two different C(/i) -corrected freely faUing observers. For any other 
pair of observers a non-zero Wigner rotation would be measured, which could lead to an 
apparent decrease in the measured EPR correlation between the two particles. 
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APPENDIX A: DERIVATION OF THE LOCAL INFINITESIMAL WIGNER RO- 
TATION EQ.(26) 

In this appendix wc derive the infinitesimal form of the Wigner rotation to 0{dT) in the 
observer's LIF given in Eq.(25) and Eq.(26), 

W\{x)^5'^^ + r~^{x)dT, (Al) 

where 

^%{x) = ^\{x)^^\{x)^0, 

A', f,r)p.(,r) -i/(,r) A-..(,x) 
■' ■' p^[x) + m 

from the formal definition of the Wigner rotation in Eq.(15) 

W\{x) ^ [L-\Kp{x)) . A{x) ■ L{p{x))] . (A3) 

As the particle moves in curved spacetime from a point from — > x" = + u°'{x)dr. 
its 4-momentum undergoes a local Lorentz transformation (LLT) p"(x) — > p'"(x) = 
A"g(x)p''(x) = p"-{x) + 5p"'{x) where 

A\{x) = + X\{x) dr, (A4) 

and 

Sp^ (x) = ^ (x) p^ (x) dr. ( A5) 

The exact form of the infinitesimal LLT A"-(x) is given in Eq.(24). However, in the subse- 
quent calculations we will only need that fact that A-g(a;) = ridcX'^i^{x) is anti-symmetric in 
its lower two indices 

Aas(^) = -hai^)- (A6) 
Note that in the observer's LIF, indices are raised and lowered with the fiat spacetime metric 
?7ac = diag(l, —1, —1, —1). The dot product of two vectors = (p°,p) and = (g°, q) is 
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given by p°-{x) qa{x) — — p{x) • q{x) where p{x) ■ q{x) = 5ijp^{x) q^{x). Further, 

because the 4-momentuin is normahzed to p°'{x) Pa{x) — m^, its time component is 
determined by it spatial 3- momentum components p'^{x) = ^/p'^{x) + m^. In the following 
we will drop the argument {x) on all spacetime dependent quantities for readability. 

In the particle's rest frame we define the standard momentum fc" = (1,0,0,0). The 
interpretation of the Wigner rotation W in Eq.(Al) is that standard boost L{p) takes the 
standard momentum /c" — > with 3-momentum components p, while the arbitrary LLT A 
takes — > p'" = A^-p'' with 3-momentum components p'. Finally, L~^{p') takes p'" — > A;", 
i.e. back to the rest momentum, which in general can differ from the standard momentum 
/c" by at most a spatial rotation k"- = H^"- k^. Thus, Wigner's little group, or the invariant 
subgroup of the massive particle's rest frame is 0(3), i.e. spatial rotations. 

The calculation of the infinitesimal form of the Wigner rotation will proceed by expanding 
Eq.(A3) to 0{dT) using the infinitesimal forms of each of the component matrices. The 
infinitesimal form of the arbitrary LLT A is given by Eq.(A4), so our first goal is to find the 
infinitesimal form of the standard boosts L{p) and L~^{p'). The form of the standard boost 
is given by Eq.(13) (with p^ — >• p"(a;)) which we repeat below 







p° 

m 








= I 
m' 




= -P± 
m' 








-(7-1) 


|.|^ = (1,2,3), 


(A7) 



where 7 = p°/m — E/m = e is the particle's energy per unit rest mass. Note that for the 
flat spacetime metric ry^g =diag(l, — 1,— 1,— l),P5=p° and pi — — p* which then agrees with 
the expressions in [17]. 

To find L~^(p'), note that L~^{p) = L(—p), i.e. the inverse of the standard boost L{p) 
which takes A;" p" is simply the standard boost in the anti-parallel direction along —p. 
However, under the inversion p — > — p, only the and - terms in Eq.(A7) which are 
linear in p change sign, while ^ and L'-. do not. As such, let us first consider L(p'), and 
then develop L~^{p') by changing the sign of the argument of the former. 
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With the substitution — > p'" — p"' -\- Sp"' in the first two fines of Eq.(A7) we fiave 

" mm 

i',(/) = ^ + < i>-) = -^-^. (A8) 
" mm' mm 

witfi given by Eq.(A5). To fiandle ) we note tfiat to 0{5p) (i.e O(c^t)) we fiave 

|p ' P = |p + 5p| ~ + 2p • 5p so tfiat 

1 1 /^^ 2p-(5p 



Tfierefore 



\ m m / |p p \ |p p 

Tfius, separating out tfie first term in eacfi of tfie equations Eq.(A8) and Eq.(A9) as com- 
posing L{p), and cfianging p — p in Eq.(A7) to form L~^{p) — L{—p) we fiave 

L-^{p') = L-\p) + M{p} dr, (AlO) 

wfiere we fiave defined tfie matrix M(p) by 

Sp^ 



M\{p') 



m 



M\{p')^-^A M\{p) = ^A 
^ m ^ m 

MK(p') = - 1) {^^^;^ - g| - ^ ^ (All) 



m 



and we fiave defined 



Sp^ ^ Sp^dT = A«.p^ = -As,. (A12) 



We can now write out tfie Wigner rotation using Eq.(A4), Eq.(A7) and Eq.(AlO) as 

W = L-\p')KL{p) 

^ {L-\p) + M{p) dr) (/ + Xdr) L{p) 

= 7 + {L-\p) A L{p) + M{p) L{p)) dr + 0{dr^) 

= I + ^dr + 0{dT^), (A13) 
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where the infinitesimal Wigner rotation is given by 

r^^=[L-\p)XL{p) + M{p)L{p)Y^^, W\ = 6'^^ + rf^dT. (A14) 

We turn next to evaluating Eq.(Al4) for the three cases (i) i}^ ^, (ii) 1?*^ and (iii) t?*-.. 
Consider the evaluation of ^ 

= L-\p)\X'^^L{p)\ + M{p)\L{p)%, (A15) 

which consists of two terms. When expanded, the first term in Eq.(Al5) becomes 

L-\p)\X'f^m^^ = L-\p)\X^.L{p)^^ + L-\p)^.X\L{p)\ + L-\p)^.y 

= 7A°.- + ^A^.7 + ^A%-. (A16) 
' m m " m ^ m 

The last term in Eq.(16) A-jp^/m^, vanishes due to the anti-symmetry of X^j. The second 
term in Eq.(16) can put in the form of the negative of the first term with the manipulations 

Pi A'o = P' Aj6 = -P' Ki = 

and therefore, Eq.(Al6) vanishes. 

The second term of Eq.(Al5) can be expanded to 

Mip)\m\ = M{p)\L{p)\ + M{p)^.L{p)\, 

= ^7+^i (A17) 



m mm 



Using 



Spi = XfPa^XfP6 + ^h^, (A18) 



Eq.(Al7) takes the form 

A°j-7 + -[Vp6 + A/p5-]^. (A19) 

Again the third term vanishes due to the anti-symmetry of A^--. Using Pq = 7717 and X~^ = 
— A°j, the second term cancels the first term, and Eq.(Al9) vanishes as well. The results of 
Eq.(Al6) and Eq.(Al7) then shows that = as stated in Eq.(A2). 

The evaluation of and proceed in an analogous fashion, albeit with considerably 
more algebra which, though lengthy, is straightforward. The following relationships prove 
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useful in manipulating these expressions into the form given in Eq.(A2). Note that p-Sp = 
S'fjP^Sp' — —pi5p^, where in the last line we have used the fact that = — p- in the flat 
spacetime metric Expanding the normalization of p'", namely (p'")^ = |p 'P + fn^ 
with p'" = p" + 5p" one easily derives the useful relationship p^5p^ = p ■ 5p. Finally, the 
normalization of the unprimed momentum (p°)^ = |p P + rn^ can be put into the form of 
a useful identity m^(7^ — l)/|pp = 1 using p° = Twy. These relations, along with the 
anti-symmetry of Aj-- and p° = pg and p* = — pj reveal after some algebraic effort that 
T^K — 1^^. — and that takes the non-zero anti-symmetric form of a rotation matrix 
given in Eq.(A2) (Eq.(26) in the text). 

APPENDIX B: DERIVATION OF THE 0{n) CORRECTED 4- ACCELERATION 
EQ.(70) 

We wish to derive the quantum mechanical correction to the 4-acceleration a°'{x) in 
Eq.(70) from the corresponding correction to the 4- velocity v°'{x) given in Eq.(67). Since 
we are only interested in keeping terms to 0{h) let us write 

v^{x) = Ua{x) + ^Sva{x) + 0{h'^), a°(x) = hSa'^{x) + 0{h^), (Bl) 

Zi I fill 

where 

5v^{x) = - {Dj^^\x)) ^p'^^\x). (B2) 

From Eq.(69) we have to 0{fi) 

= u^{x)2Dip Ua] {x) + — [ Sv'^{x)Dip Ua] (x) + uf^ix) dVa] (x) ] (B3) 

Consider the expression 2D[isUa]{x) which appears in both the first and second terms in 
Eq.(B3) 

2D[fsUa](x) = 2V[f3U^]{x) = [df),Ua{x)] = -[df),daS{x)] = 0. (B4) 

In the first equality we have used the property that the action of the total covariant derivative 
Df^ on a world vector is just the Riemann covariant derivative Vp, while in the second we 
have used the symmetry of the Christoffel symbols in their lower two covariant indices. 
Finally, we have used Eq.(53) which states that the 4- velocity is the normal to surfaces of 
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constant action Ua{x) — —daS{x). Alternatively, note that the first term in Eq.(B3) is 

U^{x)2Dli3Ua]{x) = u'^{x)Vf3Ua{x) - u'^ {x)V aUp{x) = 0, 

where the first term vanishes since u°'{x) is the tangent to the classical geodesic and the 
second term vanishes since the normalization of the classical 4- velocity u°'{x)ua{x) — 1 



The second line of the last equahty in Eq.(B5) vanishes when we invoke the equation of 
motion u^{x)D/3xjjQ^\x) which states that ijjlf\x) is parallel transported along the u^{x) 
congruence. In the first line of the second equality in Eq.(B5) we utilize the expression for 

the commutator of the total covariant derivative acting on a spinor, which is proportional 
to the Riemann curvature tensor Eq.(44) and its adjoint 



since a'^^{x) is a Hermetian matrix. Substituting these last expressions into Eq.(B5) leads 
to our final expression for the 4- acceleration 



implies that Va{u^{x)uj3{x)) — u^{x)'VaUp{x) — 0. 

The remaining third term in Eq.(B3) can be expanded as 

aa{x) = -^u'^{x) Dip6va]{x), 



= ^^uf'ix) [4''\x)Dy,D^^4-\x) - {D^^D^^4-\x))4-\x)] 

+ —uP{x) \Dpi^^^\x)D^^l^^^\x) - {Djt\^))Dpi^^^\x)\ . 
mi I ^ ' i 



(B5) 



[D^,D,]^{x) = '-R^^^s{n^))<T''\x)^P{x) 
[D^,D,]ij{x) = -'-R^,^s{V{x))ij{x)a^\x) 



(B6) 



(B7) 



aa{x) 



n_ 

4m 

h 

4m 



(B8) 



(B9) 



usmg (J 




e.^(a;) e~^{x) a^'^ and Rap^s{x) = e^^{x) e/(a;) R^pd^)- 
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